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Introduction 

The purpose of this paper and several which will follow it is to apply 
semiclassical methods developed in the Book of the author [17] to the 
theory of heavy atoms and molecules. Because of this we combine our 
semiclassical methods with the traditional methods of that theory, mainly 
function-analytic. 

In this paper we consider the case without magnetic field. Next papers 
will be devoted to the cases of the self-generated magnetic field, strong 
external magnetic field and hyperstrong external magnetic field and combined 
external and self-generated fields. Basically this paper should be considered 
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as an introduction. This paper is a revitalization of [Ivrl] but in the case 
without magnetic field. 

We explore the ground state energy, an excessive negative charge, ioniza- 
tion energy and excessive negative charge when atoms can still bind into 
molecules. 

I learned all I know about Multiparticle Quantum Theory from Michael 
Sigal, Elliott Lieb, Philip Solovej. Some remarks of Volker Bach were very 
useful for the beginner in the field. I am extremely thankful to all of them. 

. 1 Framework 

Let us consider the following operator (quantum Hamiltonian) 
(0.1.1) H = H/v:= Hv * + \ x i- x k\~ X 

1<J<N l<j<k<N 

on 

(0.1.2) fi= f\ je, ,7? = jS? 2 (R d , C) 

l<n<N 

with 

(0.1.3) H V = D 2 -V{x) 

describing N same type particles in the external field with the scalar potential 
— V (it is more convenient but contradicts notations of the previous chapters), 
and repulsing one another according to the Coulomb law. 

Here Xj e M. d and (xi, ... , x N ) e E w , potential V(x) is assumed to be 
real-valued. Except when specifically mentioned we assume that 



(0.1.4) V(x)= 



x 



l<k<M 

where > and are charges and locations of nuclei. 

Mass is equal to | and the Plank constant and a charge are equal to 1 
here. The crucial question is the quantum statistics. 

(0.1.5) We assume that the particles (electrons) are fermions. This means 
that the Hamiltonian should be considered on the Fock space Sj defined 
by (0.1.2) of the functions antisymmetric with respect to all variables 
(Xl.ft) (xa/,Qv)- 
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Here q e {1, ... , q} is a spin variable. 

Remark 0.1.1. (i) Meanwhile for bosons one should consider this operator 
on the space of symmetric functions. The results would be very different 
from what we will get here. Since our methods fail in that framework, we 
consider only fermions here. 

(ii) In this Chapter we do not have magnetic field and we can assume that 
q = 1; for q > 1 no modifications of our arguments is required and results 
are the same albeit with different numerical coefficients. In the next chapters 
we introduce magnetic field (external or self-generated) we will be interested 
in d = 3, q = 2 and 

(0.1.6) H VA = ((iV - A) ■ a) 2 - V(x) 

where cr = (oi, <J 2 , cr 3 ), are Pauli matrices. 

Let us assume that 
(0.1.7) Operator H is self-adjoint on $). 

As usual we will never discuss this assumption. 



0.2 Problems to consider 

We are interested in the ground state energy E = E/v of our system i.e. in 
the lowest eigenvalue of the operator H = H w on $): 

(0.2.1) E := infSpecH on 

more precisely we are interested in the asymptotics of E/v = E(y; Z; N) as 
V is defined by (0.1.4) and N x Z := Z\ + Z 2 + ... + Zm — > oo and we are 
going to prove that 1 ) E is equal to Thomas-Fermi energy £ TF with Scott 

5 

and Dirac-Schwinger corrections and with o(Zi) error. 
Here we use notations y = (y lf ... ,y w ), Z = (Z lt ... , Z M ). 
We are also interested in the asymptotics for the ionization energy 

(0.2.2) \ N :=-E N + E N _ V 

^ Under reasonable assumption |y m — y m >\ >• Z~i for all m ^ m' . 
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It is well-known (see G. Zhislin [Zh]) that l/v > as N < Z (i.e. molecule 
can bind at least Z electrons) and we are interested in the following question: 
estimate maximal excessive negative charge 

(0.2.3) max N-Z 

N: \ N >0 

i.e. how many extra electrons can bind a molecule?. 

All these questions so far were considered in the framework of the fixed 
positions y lf ... ,y M but we can also consider 

(0.2.4) E = E N = E(y; Z; N) = E + U(y, Z) 

with 

(0.2.5) Ln ' 1 " 



l<m<m'<M 



\Ym ~ Ym' 



and 



(0.2.6) E(Z;A/)= inf E(y; Z; N) 

yi y/w - 

and replace l/v by In = —E^ + E/v-i and modify all our questions accord- 
ingly. We call these frameworks fixed nuclei model and free nuclei model 
respectively. 

In the free nuclei model we can consider two other problems: 

- Estimate from below minimal distance between nuclei i.e. 

min |y m -y m ,| 

l<m<m'<M 

for which such minimum is achieved; 

- Estimate maximal excessive positive charge 

(0.2.7) max{Z - N : E < min E(Z m ; N m )} 

N l N M- l<m<M 

N 1 +...N M =N 

for which molecule does not disintegrates into atoms 2 ). 



2 ) One can ask the same question about disintegration into smaller molecules but our 
methods are too crude to distinguish between such questions. 
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0.3 Thomas- Fermi theory 

The first approximation is the Thomas-Fermi theory. Let us introduce the 
spatial density of the particle with the state W G ij: 

(0.3.1) p(x) = pw(x) = N J |v|/(x,x 2 x N )\ 2 dx 2 ---dx N 

where | • | means a norm in C Nq and antisymmetricity of implies that it 
does not matter what variable xj is replaced by x while in the general case 
one should sum on j = 1, ... , N. Let us write the Hamiltonian, describing 
the corresponding "quantum liquid" : 

(0.3.2) S(p) = J r(p(x)) dx-J V(x)p(x) dx + ±D(p, p), 

with 

(0.3.3) D(p, p) = Jj\x- y\- 1 p(x)p(y) dxdy 

where r is the energy density of a gas of noninteracting electrons. Namely, 
(0.3.4) r(p) = sup (pw - P(w)) 

w>0 

is the Legendre transform of the pressure P(w) given by the formula 

(0.3.5) P{w) = x lM /| +1 , X! = 2(2n)- d (d + 2)- l w d q. 

The classical sense of the second and the third terms in the right-hand 
expression of (0.3.2) is clear and the density of the kinetic energy is given by 
t(p) in the semiclassical approximation (see remark 0.3.1). So, the problem 
is 

(0.3.6) Minimize functional £(p) defined by (0.3.2) under restrictions: 
(0.3.7)! 2 p > 0, J pdx< N. 

The solution if exists is unique because functional £(p) is strictly convex 
(see below). The existence and the property of this solution denoted further 
by p TF is known in the series of physically important cases. 
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Remark 0.3.1. If w is the negative potential then 



(0.3.8) 



tre(x,x,0) « P'(w) 



defines the density of all non-interacting particles with negative energies at 
point x and 



is the total energy of these particles; here means "in the semiclassical 
approximation" . 

We consider in the case of d = 3 a large (heavy) molecule with potential 
(0.1.4). It is well-known 3) that 

Proposition 0.3.2. (i) For V(x) given by (0.1.4) minimization problem 
(0.3.6) has a unique solution p = p TF ; then denote S Tf := £(p TF ); 

(ii) Equality in (0.3.7) 2 holds if and only if N < Z := Z m ; 

(Hi) Further, p TF does not depend on N as N > Z; 

(iv) Thus 



0.4 Main results sketched and plan of the 
chapter 

In the first half of the Chapter we derive asymptotics for ground state energy 
and justify Thomas-Fermi theory. 

First of all, in Section 1 we reduce the calculation of E to calculation 
of N 1 (H M/ — v) and to estimate for D(e(x, x, u) — p, e(x, x, u) — p) where 
Ni(Hw — v) = Tr(/-/n/ — v)~ is the sum of the negative eigenvalues of operator 
Hw — v Hw = D 2 — W, 1/1/ = W TF , p = p TF are Thomas-Fermi potential and 

3 ) E. Lieb, "Thomas-fermi and related theories of atoms and molecules", [L2], pp. 



(0.3.9) 




(0.3.10) 




263-301. 
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Thomas- Fermi density respectively (or their appropriate approximations), v 
is either A a/ (A/-th eigenvalue of Hw) or its appropriate approximation and 
e(x,y, u) is the Schwartz kernel of E{y) which is the spectral projector of 
H W - 

Section 2 is devoted to the systematic presentation of the Thomas-Fermi 
theory. 

Further, in Section 3 we apply our semiclassical methods and calculate 
Ni(/"/w— v) and estimate D(e(x, x, u)—p, e(x, x, v)—p) and also |A/y— v\ where 
now v is the chemical potential (which is the Thomas-Fermi approximation 
to A/y). As a result under appropriate restrictions to A/, Z and 

(0.4.1) a := min \yj - y k \ > Z 3 
we prove that 

(0.4.2) E = £ TF + Scott + Dirac + Schwinger + o(Z§) 

and 

(0.4.3) D(p„-p TF ,p„-p TF ) = o(Zl) 

where 

(0.4.4) Scott = q Z l 



72 

"m 

Km<M 



(0.4.5) Dirac = -^(367r)igi J (p TF )J dx, 

(0.4.6) Schwinger = (36vr)M j (p TF )^ dx 

and Ml is the ground state. 

Remark 0.4-1- (i) Actually we will recover even slightly better remainder 
estimate 0{Zl' 5 ) in (0.4.2) and (0.4.3) as a > Z-h +& \ 

(ii) Condition a > Z ^ bans nuclei to be so close that the repulsion energy 
between them be much larger than the total energy of all the electrons. 
Estimates in case when this condition is violated will be also proven; 

(iii) Keeping in mind that there is no binding in Thomas-Fermi theory (and 
this statement could be quantified) one gets immediately that in the free 

_ 5 5 r 

nuclei model a > Z ,_ 2i and therefore remainder estimate 0{Z*~ ) holds. 
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(iv) Due to scaling in the Thomas- Fermi theory (see proposition 2.2.1) 
£ TF ~ qlzl = q 3 (q- 1 Z) 3 7 : Scott ~ qZ 2 = q 3 (q~ 1 Z) 2 : and both Dirac and 
Schwinger are ~ gsZ? = q 3 (q~ 1 Z)i. 

In the second half of the Chapter we apply estimate (0.4.3) to investi- 
gate negatively and positively charged systems. In Section 4 we consider 
negatively charged systems and derive an upper estimate for the excessive 
negative charge (A/ — Z) such that I a/ > and ionization energy I a/ itself. 

In Section 5 we derive upper and lower estimates for I n + v and an upper 
estimate for the excessive positive charge for which in the framework of the 
free nuclei model a < oo. 

1 Reduction to semiclassical theory 

To justify the heuristic formula E ~ £ TF = £(p TF ) and to find an error 
estimate let us deduce the lower and upper estimates for E. 

1.1 Lower estimate 

For the lower estimate we apply the electrostatic inequality due to E. H. Lieb: 

(i.i.i) Yl /i*/-**rWxi x N )\ 2 d Xl ---dx N > 

l<J<k<N J 

1 f ± 

-D(pv, py) - C / pl(x)dx 

with defined by (0.3.1). This inequality holds for all (not necessarily 
antisymmetric) functions W with ||v|/ H^^w) = 1. Therefore 



(1.1.2) (Hv|/,M/)> (Hv.x,V.V) + ^D(pv,pv)-C f pl{x)dx = 

l<j<N ' J 

J2 (Hw, Xj V,V) + \o(pv - p, p w - p) - ^D(p, p) - C I pl(x)dx 

1 i N ' " 

where (•, •) means the inner product in S) and Hw is one-particle Schrodinger 
operator with the potential 

(1.1.3) W=V-\x\- 1 *p, 
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where p is an arbitrary chosen real-valued non-negative function. 

The physical sense of the second term in 1/1/ is transparent: it is a 
potential created by a charge —p. Skipping the positive second term in the 
right-hand expression of (1.1.2) and believing that the last term is not very 
important for the ground state function v|/ 4 ) we see that we need to estimate 
from below the first term. 

Here assumption that W is antisymmetric is crucial. Namely, for general 
(or symmetric-does not matter) the best possible estimate is A/Ai where 
Ai is the lowest eigenvalue of H w (we always assume that there is sufficiently 
many eigenvalues under the bottom of the essential spectrum of H w ) and we 
cannot apply semiclassical theory. However, for antisymmetric v|/ situation 
is rather different. 

Namely, let Ai < A 2 < A3 < ... be negative eigenvalues of Hw (on 
Jt? = J2? 3 (IR 3 , C). Then the first term in the right-hand expression of 
(1.1.2) is bounded from below by 



where N(B), N^B) = Tr B~ are the number and the sum of all the negative 
eigenvalues of operator B respectively such that Spec ess (B) C 1R + provided 
A = A/v < 0; the latter assumption is equivalent to 



Applying the semiclassical approximation (which needs to be justified!) one 

gets 



(1.1.4) 




l<j<N 



(1.1.5) 



N(Hw) > N. 



(1.1.6) 



A) 



Afi(H w -\) + errori 



with 



(1.1.7) 




and therefore the lower estimate for the ground state energy is 




4 ) As we derive also an upper estimate for E we will get an upper estimate for this 
term as a bonus. 
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where error now includes both an estimate for J p^, dx and the semiclassical 
remainder estimate. 

Furthermore, applying a semiclassical approximation for the number 
N(H — A) of eigenvalues below A (and this number should be approximately 
N) one gets an equality 

(1.1.9) N = Af(H w 
with 

(1.1.10) Af(A-X):-- 

where erroro is a semiclassical approximation here. 

To get the best possible lower estimate one should pick up p delivering 
maximum to the functional 

(1.1.11) -J P(W{x) + v)dx + vN -^D(p.p) 

{y = A here) under assumptions (0.3.7) 12 and (1.1.3) as we skip all the 
errors. 

One can see that the optimal choice is the Thomas-Fermi potential and 
density. The above arguments are very standard in MQT with p = p TF , W = 
H/ TF 5 ) from the very beginning. 

On the other hand, let us consider the Euler-Lagrange equation for 
p = p TF under condition f pdx = N: 

(1.1.12) t'(p)-W = is (p>0), W=V-\x\~ 1 *p 
with the Lagrange factor u 6 \ Expressing p and integrating we get 

(1.1.13) N = M{H W - v) = J P'(W{x) + is) dx. 

Comparing (1.1.12) and (1.1.13) we get that with some error A ~ v. Substi- 
tuting to the first term in (1.1.6) A = v and v — W = —r' B (p) we get the 
lower estimate E > £ TF — error. 

5 ) Or some their close approximations. 

6 ) Called chemical potential and in contrast to A belonging to Thomas-Fermi theory. 



— A) + error 

= J P'(W{x) + X) dx 
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Remark 1.1.1. (i) Instead of (1.1.6) we will use a better estimate 7 ^; 

(ii) To minimize errors we will also recalculate (1.1.4) effectively replacing 
A by v. 

(1.1.14) X J= Yl (Aj-^) + ^A/> N^Hw -u) + uN. 

l<j<N l<j<N 

The advantage is that we even do not mess up with the semiclassical 
asymptotics for N(Hw — u). Further, one can replace here N by J p TF dx: 
(these quantities fail to be equal only for N > Z i.e. for v = 0). 

(iii) Recall that we assumed that A/v < i.e. (1.1.5) holds. In the opposite 
case 

(1.1.15) M{H W ) < N. 

we estimate the first term in the right-hand expression of (1.1.2) from below 
by Ni(Hi4/) i.e. we will get the same formula but with v = 0. 

1.2 Upper estimate 

To get the upper estimate one takes a test function ^(xi, ... , x/v) which is not 
a ground state here but an antisymmetrization with respect to (xi, ... ,x/v) 
of the product 0i(xi) • • • 0/v(x/v) where 0i, ... , 0a/ are orthonormal eigen- 
functions of H w corresponding to eigenvalues Ai, ... , A/v, provided A/v < 0. 
Namely this function minimizes the first term in the right-hand expression 
of (1.1.2). 

One can write 

( L2 -l) ^=Mdet(0 ( (x y )).. =1 N 

and it is called Slater determinant. Obviously, ||v|/|| = 1 and 

(1.2.2) pv(x)=tre N (x,x) 

where 

(1.2.3) e N (x,y)= £ 0y(x)0)(y) 

l<y'</V 

is the Schwartz kernel of the projector to the subspace spanned on {<f>j}i<j<N- 

7 ) With £ TF replaced by £ TF + Scott + Schwinger and with much smaller errori than 
for a simple semiclassical approximation. 
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Remark 1.2.1. If q > 2 then <f>j = (fij(x,q) and v|/(xi,^i; ... ;x/v,gv) is an anti- 
symmetrization with respect to (xi, <a.; ... ; x/v, <wv) of the product 0i(xi, £1) • • • 4>n{x-n, Sn)- 

Easy calculations show that 

(1.2.4) (HM/,v|/)= £ A ; + ^D(p„-p,p„-p) 

i</</v 

-^ D (P.P)-^ /7 |x-y|" 1 tre^ / (x,y)e A/ (x,y)c/xc/y. 

The first term in the right-hand expression again is equal to the middle 
expression in (1.1.4) which does not exceed 

(1.2.5) fii{H w -u) + uN+ \X N - u\ ■ \M(H W - v) - N\. 
Really, we need to consider (non-zero) terms which do not cancel in 

j<N \j<u 

and their absolute value does not exceed |A/v — v\ while their number does 
not exceed |N(/-/n/ — v) — N\. 

Again, discounting all the errors and considering semiclassical approxi- 
mation (including p v (x) ~ P'[W{x) + z/)) we arrive to a functional 

(1.2.6) -J P(W(x) + u)dx + uN-^D(p,p)+ 

^D{P'{W + v) - p, P'( W + v) - p) 

which needs to be minimized under assumptions (0.3.7)j_ 2 and (1.1.3). This 
functional differs from (1.1.11) which was minimized by the last term. One 
can prove that (1.2.6) minimizes as p = p TF , W = W TF and v is a chemical 
potential. So again we may pick them (or their appropriate approximations) 
up from the very beginning. 

Therefore in addition to a semiclassical errori 7 ) of the previous subsection 
we need to consider also semiclassical errors 

(1.2.7) D(tre(x,x,i/) - P\W + v), tr e(x, x, v) - P'{W + u)), 

(1.2.8) D(e(x,x, v) - e N (x,x), e(x,x, v) - e N (x,x)) 
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where e(x,y, u) is the Schwartz kernel of the spectral projector 0(r — H w ) 
of H w , 



Remark 1.2.2. (i) Recall that we assumed that \ N < i.e. (1.1.5) holds. 
In the opposite case (1.1.15) selecting appropriate 4>j{ x ) with j = N(H W ) + 
1, ... , N we with arbitrarily small error estimate the first term in the right- 
hand expression of (1.1.2) from above by Ni(Hw) as i.e. we will get the 
same formula but with v = and we also will need to estimate (1.2.7) with 



(ii) To make this case compatible with the case (1.1.5) we will need to 
estimate \u\ (and \N — Z|) under assumption (1.1.15); we will also compare 
S TF calculated for such v (or, equivalently, N as they are connected) and 
v = (and N = Z). 

(iii) Sure p TF and W TF depend on v (or N) but we will prove that for N — Z 
relatively small we can do all calculations as v = (and N = Z). 

(iv) If we are interested in the estimate for D(pv|/ — p TF ,pvu — p TF ) where 
is the ground state, we do not need to calculate a semiclassical error in 

Ni(/"/w — v); in fact we can simply stick with N^Hw — A) with \ = \ N 
under assumption (1.1.5) and A = otherwise. As a result in certain cases 
our estimate for D(pm — p TF , py — p TF ) will be better than the error in an 
approximation for E and we need the former rather than the latter for the 
results of the second half of this Chapter. Especially significant the difference 
will be when we introduce magnetic field. 

1.3 Remarks and Dirac correction 

Now almost everything is in framework of the theory we developed; the only 
missing is an estimate 



(1.2.9) 




and 



(1.2.10) 



v = 0. 



(1.3.1) 
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for a reasonable candidate v|/ to the ground state; one can find it in E. Lieb's 
Selecta 3 \ 

However if we want a more sharp asymptotics with Dirac-Schwinger 
terms, we need a remainder estimate o{Zi) or better; luckily there is im- 
proved electrostatic inequality due to Theorem 1, G. Graf and J. P. Solovej [GS] 
(see also V. Bach [Ba]). 

Theorem 1.3.1. Let N > eZ . Then for the ground state v|/ 

(1.3.2) E hf >E>Ehf-CZ!- 5 

and 

(1.3.3) E > E DS -CZM 
with some exponent 5 > where 

(1.3.4) E HF :=inff H F(n 

where in (1.3.4) W runs through Slater determinants 8 ^ and 

(1.3.5) E HF (^):= (Hv.xjV'V) + \v{p*,p*)- 

l<j<N 

\JJ \x- y\- 1 tre ] N {x,y)e N {x,y)dxdy, 

(1.3.6) E DS := Yl X J ~ \ D (P TF ' -° TF ) - K DiracP TF '5 dx, 

j:l<J<N;\j<0 

/tDirac = (27r)~ 3 gCj F; Cjf = (67t 2 / q 2 ) I is a Dirac constant. 

Here (1.3.2)-(1.3.5) are (1.15), (1.16), (1.8), (1.6) respectively and (1.3.6) 
is a combination of (1.12) and (2.2.12) of this paper 9 ). Actually we need 
only (1.3.3) and (1.3.6). 

8 ) Albeit not necessarily of eigenfunctions of Hyj. 

9 ) We do not have a coefficient \ in the definition of D(., .) but G. Graf and 
J. P. Solovej [GS] have. 
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As we are going to prove that the last terms in (1.3.5) and (1.2.4) coincide 
modulo 0(Z5- d ) we made a necessary step completely. 



2 Thomas- Fermi theory 

Thomas-Fermi theory is well-developed in the no-magnetic-field case. We 
cannot suggest any better reading than E. Lieb's Selecta 3 \ 

In the Thomas-Fermi theory N is a real nonnegative number (not neces- 
sarily an integer). 



2.1 Existence 

Let us recall that in order to get the best lower estimate (neglecting semi- 
classical errors) one needs to maximize 

(2.1.1) 4> m (W + v) := - J P{W + v)dx- i-||V(W- V)\\ 2 
given by (1.1.11) where we used equalities 

(2.1.2) D(p,p) = -(p, W - V) = j-\\V(W - V)\\ 2 , 

(2.1.3) p:=±A{W-V), 

\\.\\ means Jzf 2 -norm and W — > as |x| — > oo. 

On the other hand, to get the best possible upper estimate (neglecting 
semiclassical errors) one needs to minimize 

(2.1.4) <t>*(p>,v) := J (T{p')-Vp')dx+ l -V(p',p')-v J p'dx 
where 

(2.1.5) p':=P\W + u) 

and r(p) the Legendre transformation (0.3.4) of P. Recall that according to 
(0.3.5) 



(2.1.6),, P(w) = P'(w) = 



w. 
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and therefore 

(2.1.7) r(p) = 3 -{6n 2 q-^pl. 

Proposition 2.1.1. In our assumptions for any fixed v < 

(i) <t>*(W + v) is a strictly concave functional; 

(ii) <t>*(p) is a strictly convex functional; 

(Hi) <b*(W + v) < <t>*(p, v) for any p>0 and W, 

(iv) These extremal problems have a common solution W and p and 

(2.1.8) p= — A{W- V) = P'(W + v), 

(2.1.9) W = o(l) as |x| oo; 

(v) On the other hand, solution of (2.1.8)-(2.1.9) is the solution of the both 
extremal problems; 

(vi) Neither of these problem has a solution for u > 0; 

(vii) Function 

(2.1.10) A/» = J P\W + u)dx 

is continuous and monotone increasing at (— oo, 0] with Af(u) — > as 
v — > — oo and A/"(0) = Z; 

(viii) For v and N linked by N = M{v) solutions of the problem above 
coincide with p TF , W TF of the problem (0.3.6) and one can skip condition 
(0.3.6) 2 for N >Z and 

(2.1.11) S TF = cD(|/|/ TF + u ) + uN = 4>*{p TF , v) + uN. 

Proof. The proof of statements (i),(ii) is obvious; therefore both problems 
have unique solutions. Comparing Euler-Lagrange equations we get that 
these solutions coincide which yields (iv) and (iii). Proof of (v)-(vii) is also 
rather obvious. □ 
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Proposition 2.1.2. For arbitrary W the following estimates hold with ab- 
solute constants eo > and C : 

(2.1.12) e D(p-p TF ,p-p TF ) < 4>*{W TF + u) - 4>„(W + v) < 

C D(p - (J, p - p') 

and 

(2.1.13) e D(p' - p TF , p' - p TF ) < cD*(p, v ) - 4>*(p TF , u) < 

C D(p-p',p-p') 

withp= ±A{W- V), ft = P'(W + u). 

Proof. This proof is rather obvious as well. □ 

2.2 Properties 

Proposition 2.2.1. The solution of the Thomas-Fermi problem has follow- 
ing scaling properties 

(2.2.1) l/l/ TF (x; Z; y; N; q) = qlN* W TF (q*Z*x; A/ _1 Z; gi/V^y; 1; 1), 

(2.2.2) p TF (x; Z; y; A/; g) = N 2 q 2 p TF (qlz^x; N~ l Z; gi/V^y; 1; 1), 

(2.2.3) £ TF (Z; y; A/; q) = q\ N* S Tf '{N' 1 Z; gi/V3y; 1; 1), 

(2.2.4) u TF (Z; y; N; q) = q^ N^u TF (N~ 1 Z; q^ N^y; l ; l) 

where u JF = u is the chemical potential; recall that Z = (Z lt ... , Z M ) and 
y = (yi, ... , y/w) are arrays and parameter q also enters into Thomas-Fermi 
theory. 

Proof. Proof is trivial by scaling. □ 

As we can exclude q by scaling we do not indicate dependence on it 
anymore. 

Proposition 2.2.2. Let M = 1. Then the solution of the Thomas-Fermi 
problem has the following properties: 
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(i) W Tf (x; Z m ,y m ; N) and p TF (x; Z m ,x m ; N) are spherically symmetric (with 
respect to y m ) and are non-increasing convex functions of \x — y m |; 

(ii) As N = z m 

(2.2.5) H/ TF x min(Z m |x - y m |-\ |x - y m |- 4 ), 

(2.2.6) p TF x min(zl|x - y m | i , |x - y m |- 6 ) 

_ 1 4 

with the threshold at |x — x m | x r m = Z m 3 when W Tf x Z m and p TF x Z m . 
(Hi) As eZ m < N < Z m 

(2.2.7) -z/x|Z m -A/|s 

and (2.2.5) holds as |x — y m | < r m x — z/|Z m — A/| _1 w/iz/e 

(2.2.8) 1^/ TF x (Z m - A/)|x - y m |- 1 as \x - y m | > |Z m - A/ 1 " i ; 

fwj Meanwhile p TF = as |x — y m | > r m and p TF = 0(|Z m — A/| 2 ) as 
|x - y m | < r m . 

For atoms let ? m denote the exact radius of support of p TF . 
Proposition 2.2.3. As M > 2 

(2.2.9) ^ e H/j; F (C(x-y m ))< 14/ TF < C^H/: F (e(x-y m )) 

m m 

where W^ F stays for an atomic solution with the charge Z m located at y m 
and the same chemical potential v . 

Proof. Proof is due to comparison arguments of E. Lieb, J. P. Solovej and 
J. Yngvarsson [LSY1, LSY2]. □ 

Proposition 2.2.4. (i) Let N < Z and 

(2.2.10) £(x) = ^min|x-y m |, 

(2.2.11) cm= {w T "( x )y- 
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Then 

(2.2.12) CM<CM 



i 



ZU(x)s as £(x) < Z^a, 

C \ £(x)- 2 as Z-3 < £{x) <\Z- N\~K 

\Z-N\h(x)~* as £{x)>\Z - A/|-i; 

£(x) and £(x) ore £-admissible and 

(2.2.13) |D a l/l/ TF (x)| < C a C(x) 2 £(x)~ |a| Va : \a\ < 3, 
and 

(2.2.14) \D a W TF (x) - D a W TF (y)\ < Q((x) 2 £(x)~ 7 2 \x - y|i 

Vx, y : |x — y| < e£(x) 

^ Unless ((x) x (— estimates (2.2.14) hold for all a. 

Proof. This proof is rather obvious corollary of the Euler-Lagrange equa- 
tion. □ 

Remark 2.2.5. Let 

(2.2.15) Z m xW Vm. 
Then C(x) x C(x). 

Theorem 2.2.6. 10 ^ Consider S TF and 

(2.2.16) f TF := £ TF + L/, 

(2.2.17) ty=W(Zi Z M ; yi y M )= ]T r^T^T" 

l<m<m'<M 1 1 

Select a nucleus y m and a unit vector n such that 

(2.2.18) (y/<-y m ,n)<0 V/c 

and plug y m + an instead of x m mto £ TF and mto £ TF n ). TTien 
10 ) Theorem 1 of R. Benguria [Be] . 

n ) So all other nuclei are confined in half-space and y m moves away outside 
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(i) Sj F is a non-increasing function of a > 0; 

(ii) S^ F is a non-increasing function of a > 0; 

(iii) -(iv) For fixed a > both Sj F — £j F and £j F — 8q F are non- decreasing 
functions of N. 

Equality 

d£ TF 

(2.2.19) 

implies that (iii)-(iv) is equivalent to 

(2.2.20) u a > i/ . 
Theorem 2.2.7. 12) 

(i) For fixed Z lt ... , Z M ; y x , ... , y M and N = Z 

(2.2.21) A^^Z; Ay; N) = ^(A 3 ^; y; A 3 A/) 

is positive non- decreasing function of A > and has a finite limit as A — > 
+oo; 

(ii) This limit does not depend on Zy, ... , Zm- 

Scaling property without assumption N = Z holds for S TF and U as well. 
These two theorems imply immediately 

Proposition 2.2.8. Let (2.2.15) be fulfilled. Then 

(2.2.22) ^(Zjyi/V)- min V £ TF (Z m ; N m ) > e min(a~ 7 , Zs) 

A/i /V/vc ^— 



N 1 +...N M =N !< m < M 



(2.2.23) a = - min |y m - y f 

Z m<m' 

12 ) (1.8)-(1.9) of H. Brozis and E. Lieb [BrL] . 
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Proof. In virtue of theorem 2.2.6(i) it suffices to prove proposition for M = 2 
(all other nuclei could be pulled to infinity), and in virtue of theorem 2.2.6(iii) 
it suffices to prove proposition for Z = N. 

Then the proof is due to theorem 2.2. 7(i) and (2.2.15), which provides 
uniformity. □ 

Remark 2.2.9. In virtue of (2.2.19) the minimum (with respect to N\, ... , Nm) 
in the sum in the right hand expression is reached when vj = for all j, k. 
The same is true for a system of isolated molecules. 

Proposition 2.2.10. Let Q denote Thomas-Fermi excess energy which is 
the left-hand expression of (2.2.22). Then 

(2.2.24) D(p TF -p TF ,p TF -p TF )< CQ, p TF := £ p JF . 

l<m<M 

Proof. We follow "non-binding" proof due to Baxter (see E. Lieb S electa 3 ^). 

According to Baxter lemma there exist g, < g < p TF and h = p TF — g 
such that g * \x\~~ 1 = Vf a.e. when h > and g * |x| _1 < Vf a.e. when h = 0. 
Here V m = Z m \x - y m | _1 . 

Let a = j g dx, (3 = J h dx and let Sj F , S TF ' be Thomas-Fermi energies 
for the first atom and for the rest of molecule respectively and p TF , p TF/ be 
corresponding Thomas-Fermi densities. Then 

(2.2.25) n min^^/VO + S TF '(N')) < £,{a) + £'(/3) < 

£ 1 (g) + S\h) - eD(g - pJ F , g - p[ F ) - eD(h -p',h- p') < 
£(g + h) + J h(V 1 -g*\x\- 1 )dx- J(V 1 -\g*\x\- 1 )d\i' 

- eD(g - p[ F , g - pY) ~ eD(/7 - p TF/ , h - p TF/ ) 

where \x! are measures with densities Z 1 8(x— yi) and J2 2 <m<M Z m b{x—ym) 
and we used the superadditivity of r(p) = ps. The last expression doesn't 
exceed 

(2.2.26) £ TF - eD(g - pJ F , g - pJ F ) - eD(h - p TF/ , h - p TF/ ). 
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Using induction with respect to M we arrive to 

(2.2.27) D(p TF - pY - P TF/ , P TF - pY ~ P TR ) < 

2D(g - pY, g - pY) + 2D(/i - p TF/ , h - p TF ') < CQ 

and finally to (2.2.24). □ 

Problem 2.2.11. Find the stronger lower bound in (2.2.22) as N < Z. 
Would be the left-hand expression x min(a~ 7 + |Z — A/| 2 a _1 , Zs) ? 



3 Application of semiclassical methods 
3.1 Asymptotics of the trace 

In this subsection we calculate asymptotics of Tr(Hw — v)~ . Here we need 
to consider both inner and outer zones. 

An inner zone (near nucleus x m ) is a ball where Vj = Z m \x — y m | _1 
dominates W — V m . For a single nucleus (M = 1) it is defined by 

(3.1.1) |x-y m |<eZ m ^ 

but in the case M > 2 there are another restrictions 

(3.1.2) |x - y m | < e min [Z m (Z m + Z m /) _1 |y m - y m /| ) 

and 

(3.1.3) |x - y m | < Z m z/ _1 
but we shrink this zone to 

(3.1.4) |*-y m | < r m := e min (Z m Z" 1 a,Z m 3 ). 

Let us consider contribution of the zone X m described by (3.1.4) to 
Ni(H w — v) = Tr(H w — , both into the principal part of asymptotics and 
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the remainder. Let ip m be a partition element concentrated in X m and equal 
to 1 in {x : |x — y m | < \r* m }. Then according to Theorem 12.5.8 of [17] 

(3.1.5) Tr((Hw - v)-^ m ) = J Weyl 1 (x)^ m (x) dx + Scott m + 0{R m ) 

where Weyl^x) and Scott m are calculated for the case q = 1 and then 
multiplied by q 13 ^: 

(3.1.6) \Ney\ 1 (x):=-^- 2 (W(x) + u)l 
while R m is C( 2 ((£) = C( 3 £ calculated on its border i.e. 

(3.1.7) R m = czl + czlrj < cz\ + czh-l 

Really, one needs just to rescale x \-t (x — ym)^ 1 and r \-t rZ^r m and 

ii 

introduce a semiclassical parameter h = Z^r™- 

Remark 3.1.1. (i) These arguments work only as r m > Z m 1 (i.e. Z m > 

(ii) As > a _1 Z but a > Z~ x we define r m = a^Z - ? and we do not include 
Scott m 14 ) into the principal expression; moreover, in this case we include X m 
into a singular zone and use variational methods to estimate its contribution 
into the principal part of asymptotics; it will not exceed CZa 1 < CZ2a~*] 

(iii) Furthermore, as a < Z we set r m = Z -1 and we do not include any 
Scott m u ' into the principal part of asymptotics and include all X m into 
singular zones; using variational methods we estimate their contributions 
into the principal part of asymptotics by CZ 2 . 

Therefore 

(3.1.8) The total contribution of all inner zones into remainder does not 
exceed the right-hand expression of (3.1.7) as a > 2 ' 1 and CZ 2 as a < Z 1 . 

13 ) As operator H m — v is nothing but q copies of A — ( W + v). 

14 ) However it will be less than the remainder estimate, so we can include it into the 
principal part of asymptotics anyway. 
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Let us consider contributions of the outer zone Xq which is complimentary 
to the union of inner zones. Then 

(3.1.9) Tr({H w - z/T^o) = f Weyl 1 (x)^ (x) dx + 0{R Q ) 
with Weyl^x) defined by (3.1.6) where 

(3.1.10) Rq= I C((xfr 2 dx 
and 

_ f Ziri as £ < z-3, 

3.1.11 (=< ~ x 

\r 2 as £>z-i. 

We justify (3.1.9)— (3.1.11) a bit later by an appropriate partition of unity. 
One can see easily that contribution of {x : £{x) < Z~3} into expression 
(3.1.10) does not exceed the same expression as in (3.1.8) and that contribu- 

1 - 5 

tion of {x : £{x) > Z~$} into (3.1.10) does not exceed CZs. Then we arrive 
to 

Theorem 3.1.2. Let W = W TF , N^Z, W = W TF . Then 
(3.1.12) Tr{H w -i/)"= / Weyl^x) dx + Scott m + O(R), 

J Km<M 



with 
(3.1.13) 



CZ3 + CZza 2 as a > Z 1 , 

CZ 2 as a< Z- 1 . 



Proof, (i) Consider v = first. Then we just apply ^-admissible partition 
of unity. Sure, (£ < 1 as £ > 1 but we can deal with it either by taking 
( = £~ x here or considering it as a singular zone and applying here variational 
estimate as well. 

(ii) Variational estimate works for v < as well; furthermore, zone {x : 
W(x) < (1 — e)v, £ > is classically forbidden. 

(iii) However as v < — c we have a little problem as W TF is not very smooth, 
it is only ^2 as W x —v. The best way to deal with it is to take 

15) F ronl the point of view of generalization to the case when magnetic field is present 
and it is not too weak, so magnetic version of W TF has multiple singularities. 
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e£-mollification with e = h 1 s , K = ((£) 1 , use rough microlocal analysis 
of Section 4.5 of [17] and framing; it will bring an approximation error not 
exceeding Ce^hr' i \v\ which does not exceed \v\ = 0(\Z — A/ 1 3 ) . We leave 
easy details to the reader. □ 

Now we arrive to the lower estimate for E: 
Corollary 3.1.3. isWxZ 

(3.1.14) E > £ TF + Scott — CR 

with R defined by (3.1.13). 

Proof. We know from subsection 1.1 that 

1 5 

(3.1.15) E > Ni(Hw -v) + uN- -D(p, p) - CZ* 
as W is given by (1.1.3). In virtue of theorem 3.1.2 

(3.1.16) E > J Weyl^x) dx + u J Weyl(x) dx - ^D(p, p) + Scott - CR 

where we also plugged instead of N as N < Z (and v < 0) 

(3.1.17) N = J \Ney\{x)dx 
with 

(3.1.18) Weyl(x) = -^{W{x)+v)\. 

One can check easily that three first terms in the right-hand expression of 

(3.1.18) constitute exactly <b*(W + v) coinciding with £ TF as 1/1/ = 1/1/ TF .D 

Remark 3.1.4- (i) As a < Z ' J using the same method one can prove a 

3 1 

slightly better remainder estimate-with Z^a~^ replaced by 

(3.1.19) J2 ™n({Z m + Z m ,) 3 2\x m -x m ,\-K(Z m + Z m ,) 2 ^) 

l<m<m'<M 

allowing to lighter nuclei to be closer one to another. 

(ii) To improve it further allowing lighter nuclei to be closer to heavier ones 
one needs to improve Theorem 12.5.8 of [17] which seems to be too difficult 
task for a such little gain. 
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3.2 Upper estimate for E 

Recall that in virtue of subsection 1.2 

(3.2.1) E < Hi(H w — v) + isN — ^D(p, p)+ 

\X N - v\ ■ |N(Hw - v) - N\ + ^D(tre/v(x,x) - p, tr e N (x, x) - p) 
and thus we need to estimate two last terms in the right-hand expression. 

3.2.1 Estimating |A/v — v\ 

First we need to estimate |A/v — v\. We will use the heuristic equality 
N(Hw — A/v)~A/or more precisely 

(3.2.2) ! 2 M{H W - X N - 0) < N < H{H W - X N + 0) 
and equality 

(3.2.3) J Weyl(x) c/x = min(/V, Z) 
where the (3.2.2) 2 is valid only if \ N < i.e. N(Hn/) > N. 
Case Aw < v. 

Then we will use (3.2.2) 2 and to calculate U(Hw — A/v + 0) we will use 
semiclassical approximation: 

(3.2.4) M{H W - A/v + 0) = J Weyl(x, \ N ) dx + O(/? ) 
with the semiclassical error 

(3.2.5) Rq = J fr 1 dx 
with integral not exceeding 



7 , 

^{|x|>Z-3} 



CZ / |x|- 2 c/x+C/ , |xr 5 c/xxCZf; 

J{\x\<Z-\} 
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again we need to consider separately cases of N > Z and v = when integral 

(3.2.5) is taken over IR 3 and N < Z and v < when it should be taken 
over {£{x) < C(Z - N)- 3}; in the latter case to cover non-smoothness we 
consider an approximation (mollification) of W and an approximation error 
£2 fir 3 <^ 1. Therefore 

(3.2.6) J Weyl(x, \ N ) dx > N - CZ^ 
comparing with (3.2.3) we conclude that 

(3.2.7) J ((H/(x) + {W{x) + X N )lj dx < CZl 

Here an integrand is non- negative (since A a/ < v) and one can see easily that 
the main contribution to integral is delivered by zone {x : £(x) x |Aa/|~ 5 } 16 - ) 
and the whole integral is x | A/v | ^ | — \n\- Therefore (3.2.7) yields 

|A/v - is\ < C(\u\ + |A/v - z/|)?Z5 
which is equivalent to 

(3.2.8) \X N -v\< CZl + C|i/]3Zi x CZ§ + C(Z - N)\zs < CZ. 
In particular 

(3.2.9) If H < c Z§ (i.e. (Z - N)+ < c x Z%) then |A W | < C Z§ 
and 

(3.2.10) |A/v - i/| • |N(H W - v) - N\ < CZ ■ Z§ = CZL 
Case A/v > ^. 

Then we will use (3.2.2) 1 but if N(Hw) < A/ then integral (3.2.5) is diverging. 

To avoid all related difficulties we will consider first the case when 
we necessarily conclude that |A/v| > (1 — e)M- To do so note that even 

!6) Provided |A/v| < QZs . On the other hand, if |A«| > QZs one can see easily that 
(3.2.7) is x Z because \v\ < cqZs due to N x Z. 
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if the main contribution to integral (3.2.7) is delivered by zone {x : 
£ x (Z — A/^A/vl" 1 x \v\ * | A a/ | 1 j" , we will ignore it considering a zone 
{x : £(x) < Cb|H"*} instead of {x : £(x) < C |z/| ? | A a/ | ~ 1 > 18) - 
One can prove easily that 

(3.2.11) The contribution of zone {£ < Co\v\~*} to the semiclassical remain- 

2 

der when calculating N(Hw — A/v) does not exceed CZ3. 
Therefore we arrive to 

(3.2.12) \\ N -u\< C\v\*zl x C(Z - N)\Z* < CZ 
(cf. (3.2.8)). In particular 

(3.2.13) If \v\ > CZl (i.e. (Z - N) > CZ§) then \X N \ x \u\ 

(cf. (3.2.9)). Then one can easily recover (3.2.8) completely. Since \N(Hw — 
v)-N\< CZ\ we arrive to (3.2.10). 

Case \v\ < 7] = CZv. 

3 2 

This case (i.e. (Z — N) < r\* = CZ?) is the most important one. The easiest 
way to tackle it is to pick up p TF and W TF calculated as if v = i.e. Z = N; 
that means the change of the test function W in the upper estimate. 

We need to modify an upper estimate of Xa<y</v -V To do this we note 
that 

(3.2.14) A/v > -77 
and 

(3.2.15) A number of eigenvalues between \ N (if \ N < 0) and does not 

2 3 

exceed CZs + C??* 

which can be proven easily by our standard methods. Therefore 

(3.2.16) J2 X J < JfH w + C v{Z~i + V 1 ) 

l<j<N 

5 20 

and the last term is less than CZ3 as long as r] < CZ 21 which is fulfilled. 
In this last case 19 ) we arrive to an upper estimate with £ TF calculated 

17 ) Now an integrand is non-positive. 

18 ) Obviously these zones coincides as \u\ x |A/v|. 

19) Pending an analysis of the next subsubsection. 
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as if v = i.e. £ TF (Z; y; Z) which is less than £ TF (Z; y; A/). Actually the 
difference between these two is x < 1 77 j 5 . 

3.2.2 Estimating D-term 

We need to estimate the last term in the right-hand expression of (3.2.1); 
we estimate it by 

(3.2.17) C D(tre( x, x, v) — P'(W + u), tr e(x, x, z/) — P'(W + z/)) 
+ C D(tr e(x, x, z/) — tr e w (x, x), tr e(x, x, z/) — tr e w (x, x)) 

+ C D(p- P'{W + u),p- P'(W + u)) 

where the last term vanishes as p = p TF , W = 1/1/ TF ; however for technical 
reasons we want to avoid this assumption. 

Estimating the first term. 

To estimate the first term in (3.2.17) we apply semiclassical asymptotics 

(3.2.18) tr e(x, x, u) = Weyl(x) + 0(( 2 r 1 ) 

where Weyl(x) = P'(W(x) + v) and therefore this term does not exceed 

(3.2.19) J j C(x) 2 C(y) 2 ^(x)- 1 £(y)- 1 |x-y|- 1 c/xc/y. 

Estimating this integral by the double sum of integrals over domains {(x, y) : 
£(x) < Z-lj(y) < Z-\} and {(x,y) : £(x) > Z~\,i{y) > Z~\} we get 

(3.2.20) CZ 2 ff i \x\- 2 \y\~ 2 \x - y]- 1 dxdy 
and 

(3.2.21) C [f , \x\- 3 \y\- 3 \x - yl' 1 dxdy 

J J {\x\>Z~l ,\y\>Z-\} 

respectively and rescaling we get the same integrals but both with the 
"threshold" 1 rather than Z~5 and both with factor Zs rather than Z 2 or 1 
respectively; one can see easily that both integrals (without factor Zs) are 
x 1 and then (3.2.19) is 0(Z§). 
Therefore we proved 
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Proposition 3.2.1. As W = W TF the first term in (3.2.17) does not exceed 

czl 

Remark 3.2.2. (i) For N > Z and u = we used that C < Q^~ 2 for 
£> Z-i; 

(ii) for N < Z and v < we used that zone {x : £ > C(Z — A/)~3} is 
classically forbidden ( W + v < there) and therefore integral is taken over 
zone {x : i < C(Z - A/)-|} where ( < Qi' 2 ; 

in As N < Z W TF is not very smooth near W + v = but one can handle 
it by rescaling arguments; 

(iv) Alternatively (preferably 15 -*) one can replace W TF by its mollification 
WJ F . 

Estimating the second term. 

Consider now the second term in (3.2.17). Due to arguments of the previous 

5 

paragraph modulo O(Zs) one can rewrite it as 

(3.2.22) C 1 D(P'(W + X N ) - P'(W + u),P'(W + \ N ) - P'{W + u)). 
Really, if we replace tr e(x, x, v) and tr e/v(x, x) by Weyl(x, v) and Weyl(x, \ N ) 

5 

respectively we make a semiclassical errors estimated by Zs provided either 
A/v<^orA/yXz/ which is always the case unless \u\ < CZ$ but in this case 
we just "cheat" resetting everything to the case v = 0. 

Let us estimate (3.2.22). According to the previous subsubsection there 
are two cases: 

(i) N > Z - CZl, in which case \u\ < CZl and |A W | < CZ§, and (3.2.22) 
does not exceed 

CW-\n\ 2 II \x-y\- 1 C(x)C(y)dxdy 

J J {e(x)<L,£(y)<L} 

+C ff \x - yrWCCr) 3 dxdy 

J J {e(x)>L,e(y)>L} 

with L = | A/v — v\~^\ one can calculate easily that both terms are x 
C A/v - v U < CZt < Zs. 
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2 

(ii) N < Z — CZ5, in which case (3.2.22) does not exceed 

CW-Xn\ 2 [[ \x-y\- 1 C(x)C(y)dxdy 

J J {t(x)<L,i{y)<L} 

with L = \v\~* and this integral does not exceed |A/v — v\ 2 \v\~^ which due 
to (3.2.8) does not exceed C(Z - N)lzl < CZ\. 

Estimating the third term. 

The third term in (3.2.17) does not exceed 

C 1 D(p-p TF ,p-p TF ) + 

^(P'iW + u) - P'(W TF + u),P'(W + u) - P'(W TF + 

and we leave to the reader an easy proof that both terms here are 0(Z^~ S ) 
as W is a described mollification of W TF . 

3.2.3 Finally, a theorem 

As we finished an upper estimate for E we arrive to 

Theorem 3.2.3. Let N x Z and let W be a ground state. Then 



Z3 + Ca 2Z2 as 3> Z \ 

Z 2 as 3< Z 1 



(3.2.23) |E-£ TF -Scott| < C 

where a is the minimal distance between nuclei. 



3.3 Improved asymptotics 

15 

So far as a > Z 3 we recovered only 0(Z5) for both error estimate in E and 
(as a coproduct, see subsection 3.4) for D(pvu — p Tf , pm — p TF )- 

5 _1 

Our purpose is to improve them to o(Zs) (or slightly better) as a > Z 3 
and recover Schwinger and Dirac terms. To do so in the lower estimate for 
E one just need an improved electrostatic inequality (see Theorem 1.3.1) 
and also improved semiclassical estimates in Tr(H w — v)~ and U(H W — v). 
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For improved upper estimate we will need also to improve estimate 
|D(pm/ — p TF , D(pm/ — p TF ) for the test function W and apply an estimate 

(3.3.1) \±JJ \x-y\- l tre ] N {x,y)e ] N (x,y)dxdy+ 

which is due to Theorem 6.3. 1? 20 -* of [17]. 
Remark 3.3.1. (i) Only contribution of zones 

(3.3.2) {x : Z~Hi < | x _ y m | < Z~h Sl } 

with b := min(aZ5, Z) should be considered as contribution of both zones 
{x : £(x) < Z-lb- 5 *} and {x : t{x) < Z^i} are 0{Zlb- 5 ). 

(ii) Only m with Z m > Zb~ &2 should be considered as contribution of others 
is lesser than the remainder estimate. 

Proposition 3.3.2. Let $ be a (small) parameter such that b~ &2 < i? < 1; 
consider m with Z m > Z$ S K Let<p{x) = ip{r^(x-y m )) with if) E ^{B(0, 1)) 
and r defined below. 

(i) Further, let (Z — N)~s > r = Z~^"&~ x . Then inequalities 

(3.3.3) \[<f>{x)f (e(x, x, A') - Weyl(x, A')) dxd\' 

J J — OO 

I 5 A 

— Scott — Schwinger] < CZ^§ , 

(3.3.4) | J 0(x)(e(x,x, A) - Weyl(x, A)) dx\ < CZU 5 
and 

(3.3.5) D(0(e(x,x,A) - Weyl(x, A)) , 0(e(x, x, A) - Weyl(x.A))) < 

czU s 

hold with some exponent 5 > for all A < and for <fi which is r-admissible. 

20 ) This theorem implies the above estimate with 5 — 1 which is definitely overkill for 
our purposes. 
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(ii) On the other hand, let (Z — N) 3 < r. Let W be a constructed above 
mollification of W Tf . Then 

(a) Estimates (3.3.3) -(3.3.5) hold for all A < v; 

(b) Further, estimates (3.3.4), (3.3.5) hold for A G [za 0] such that N(H — 
A) < N. 

(c) Furthermore, in this last case 

(3.3.6) \\-v\<CZl(Z-N)-*$ 5 . 

To prove these statements we need to study behavior of the Hamiltonian 
trajectories. First we want to prove that in the indicated zone 1/1/ TF is a 
weak perturbation of Wjf which is a single atom Thomas-Fermi potential 
with Z m and with z/ m = v. 

Proposition 3.3.3. In the framework of proposition 3.3.2 in B(y m , r) 

(3.3.7) \D a (W TF - W J J)\ < c a W J J\x - y m |- |c V. 

This estimate holds for all a as l/l/ TF /(— u) is disjoint from 1; otherwise it 
holds for I a I < 3 and 

(3.3.8) \D a {W TF - Wl ¥ )(x) - D a {W TF - W^){y)\ < 

cW^\x -x m \-i\x - y\h 5 

as \a\ = 3 and \x — x m | x |y — y m \ x (Z — A/)~s . 

Proof. An easy proof based on variational approach is left to the reader. □ 

Next let us consider manifold Z A = {(x, £) : H(x, £) = A} with H(x, £) = 
|£| 2 — H/(x), and let us introduce a measure ]i\ with the density dxd^ : dH 
on Za; this measure is invariant with respect to the Hamiltonian flow with 
the Hamiltonian /-/(x,£). Note that ^(^a) x Z. 

Proposition 3.3.4. In the framework of proposition 3.3.2 there exists a 
set T! x C T\ such that 

(3.3.9) jt A (r' Ai ,) < C$ S Z 
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i 2 

and through each point (x, £) belonging to i?(Z~3, Zs) -vicinity of T\ \ Z' A l? 
£/iere passes a Hamiltonian trajectory (x(t),£(t)) 0/ H of the length T = 

(3.3.10) I D a 1 2 (x(t)Z3,^(t)Z-3)| < GT K Va:|a|<m 
and 

(3.3.11) |x(t) - x(0)|Zs + |f (t) - f (0)|Z-§ > $ K \t\Z 
where m is arbitrary and K, 5 depend on m. 

Proof. We will just sketch it. First, let us consider the case M = 1. In 
this case we first include into E Al? all the points (x, £) with trajectories not 
residing over {$ 5 'Z"3 < £(x) < Z~5i?~ 5 '} for all t : \ t\ < 7~ 21) . 

Further, if (Z — A/) + > Z& we also include into T! x $ all the points (x, £) 
with trajectories not residing over e(xi,(l-#')r); here r is an exact radius 
of suppp TF . Then (3.3.9)-(3.3.10) hold and we reduce 5, 5' if necessary. 

It is known that not all the Hamiltonian trajectories are closed (see 
Appendix 6.B). Then one can prove easily that adding to Z Al? the set 
satisfying (3.3.9) we can fulfill (3.3.11) as well. One can find the similar 
arguments in the proof of Theorem 7.4.12 of [17]. 

The general case is due of this particular one, proposition 3.3.3 and 
trivial perturbation arguments. □ 

Proof of proposition 3.3.2. Now statements (i), (ii) follow from proposi- 
tions 3.3.3 and 3.3.4 and the standard arguments; as (Z — A/) + > ZD 5 we 
need to mollify 1/1/ TF in the standard way. 

To prove statement (iii) we can use decomposition (16.3.1) of [17] and 
apply to the contribution of zone {(x,y) : |x-y| > Z~ *■&*'} the same 
standard arguments. Meanwhile one can notice that contribution of zone 
{(x, y) : |x - y| > Z^ 5 '} is C^ZM). □ 

Combining all these improvements we arrive to 

Theorem 3.3.5. Let N x Z. Let a > Z "3 and let W be a ground state. 
Then 

(3.3.12) |E -£ Tf - Scott - Dirac - Schwinger| < CZ§ (Z~ 5 + {aZ^y 5 ) . 
21 ) I.e. trajectories, leaving this "comfort zone" for some t : |t| < T. 
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3.4 Corollaries 

3.4.1 Estimates for D(py — p Tf , py — p Tf ) 

Recall that in the lower estimate there was in the left-hand expression a 
non-negative term ^D(p^ — p TF , p-y — p TF ) which we so far just dropped. Also 

1 3 

recall that the term a^zZs in the remainder estimate appeared only because 
we replaced Tr(H w — v\~ by its Weyl approximation (with correction terms) 
which we by no means need for estimating this term as Tr(H w — v)~ was 
present in both lower and upper estimates. Then we arrive to 

Theorem 3.4.1. Let N x Z and let ^ be a ground state. Then 

(3.4.1) D(p„-p TF ,pv-p TF ) < CQ := 

{ 5 1 
Zs as a < Z 

zl(Z- 5 + (aZ5)- 5 ) as a>Z~\. 

3.4.2 Estimates for distance between nuclei in the 
free nuclei model 

Let us estimate from below the distance between nuclei in the stable molecule 
in the free nuclei model (with the full energy optimized with respect to the 
position of the nuclei). 

Theorem 3.4.2. Let M > 2 and let Z m ^Z for all m be fulfilled. Assume 
that 

(3.4.2) E(Z;y;A/)+ ^ ZjZ k |y m - y^" 1 < 

l<m<m'<M 

min H z ^N m ) 

Nl N M- l< m <M 

N! + ...+N M =N 

Then 

(3.4.3) |y m -y m '| > z~ k+Sl y m ^m' 

and 

(3.4.4) |E TF (Z; N) - ^ TF (Z; N) - Scott - Dirac - Schwinger| < CZH : . 
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Proof. Note first that |E| < CZs and in virtue of theorem 3.2.3 we can 
replace E by £ TF with 0(Z 2 ) error: 

(3.4.5) £ TF (Z;y;A/) + ^ ZjZ k |y m - y^" 1 < 

l<m<m'<M 

min ^ £ TF (Z m ;A/ m )-CZ 2 ; 

Wl N M- l< m <M 

N! + ... + N M =N 

which in virtue of proposition 2.2.8 is impossible unless a~ 7 < CZ 2 i.e. 
a > eZ~ 2 7. 

Then, again in virtue of theorem 3.2.3 we can replace E by S TF + Scott 
with 0(Z5) error. Let us take into account that for molecule Scott equals the 
sum of Scott m . Therefore in (3.4.5) we can replace CZ 2 by CZs. Applying 
again proposition 2.2.8 we conclude that a > eZ~2i. 

Let us improve this estimate further. In virtue of theorem 3.3.5 again 
we can replace E by £ TF + Scott + Dirac + Schwinger with 0(Z3~ S2 ) error. 
However one needs to compare Dirac-Schwinger correction for molecule with 
the sums of such corrections for separate atoms: 

Proposition 3.4.3. If a > Z s +<5 i and 

(3.4.6) £ TF (Z;y;A/) = £ £ TF (Z m ; N m ) + 0(z5"*) 

l<m<M 

where N = Ni + ... + N M then 

(3.4.7) f(p Tf ) l dx= J2 [(pl^dx+OiZ^) 

where p T ^ = p TF (x — y m ; Z m ; N m ) are atomic Thomas- Fermi densities. 

Proof of proposition 3.4.3 will be provided immediately. Therefore Dirac 
and Schwinger correction terms for molecule are equal to the sums of Dirac m 
and Schwinger m correction terms respectively with 0(Zi~ &2 ) error and in 
(3.4.5) we can replace CZ 2 by CZs - * 52 . Applying proposition 2.2.8 again we 
conclude that a > Z^ b . □ 
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Proof of proposition 3.4.3. Note that the left-hand expression of (2.2.24) is 
equal to 

(3.4.8) ||V(l/1/ TF - W TF )\\ 2 , W TF : = ^ W^ F 

l<m<M 

and therefore this expression is less than CQ < CZ^~ Sl as well. Then as a > 
Z~i + 1 we conclude that if we restrict norm to zone {x : |x — y m | 
we can replace p TF and W TF by p]J and Wjf respectively in (2.2.24), (3.4.8). 
Using Thomas-Fermi equations we conclude that in this zone 

(3.4.9) \D a {W TF - W^)\ < C a Wl F t^Z- 52 Va : \a\ < 2 

and then |(p TF )3 — (Pm F )^| < C(pl^)^ Z~ Si which implies (3.4.7) because 
f(Pm F ) 5 dx x Zs and contributions of zone {x : £(x) > Z-3+ 5 } to each 
integral is 0(Z^~ 5 ). □ 



4 Negatively charged systems 

In this section we consider the case N > Z and provide upper estimates for 
excessive negative charge (A/ — Z) as I a/ > and for ionization energy \ N . 

4.1 Estimates of the correlation function 

First of all we provide some estimates which will be used for both negatively 
and positively charged systems. Let us consider the ground-state function 
V(xi, qi, ... ; x/v, gv) and the corresponding density py/(x). 

The crucial role plays estimate (3.4.1) D(p^, — p 1F , — p TF ) < CQ of 
theorem 3.4.1 and the difference between upper and lower bounds for E/v 
(with N 1 (/-/ M/ — v) + uN not replaced by its semiclassical approximation). 

Let us consider the classical density of the electron system 

(4.1.1) g x _(x) = b ( x - x j) 

1<J<N 

and the smeared classical density 

(4.1.2) Qx,e{ X ) = Qx*4>e= ^2 ^ £ ( x ~~ X J') 

l<j<N 
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where e will be chosen later; here (x, q) = (xi, si; ... , x N , s N ) G (IR 3 x C <7 ) /v , 

(4.1.3) e (z) = £ _3 0(z£ _1 ), G < ^ oo (6(0, |)) is a spherically symmetric, 
non- negative function such that f 0(x) dx = 1. 

Then / e (x)f(x) c/x = f(0) + 0(e 2 ) for any f G <^ 2 . Let us consider 

(4.1.4) K N {x) := l D (fo(.) -p(-), &(•) 

where p = ^A(W— V), W is either W TF if z/ = or a "good" approximation 
for it, constructed in the previous section. 

Using Newton's screening theorem 22 ) we conclude that 

(4.1.5) J2 \ X J ~ x ^ * ^D( fe (-), &(•)) - C ^ lA/ 

l<J<k<N 

where the last term estimates 

D(<f> e (x-Xj),<f> e (x-Xj)). 

l<j<N 

On the other hand, 

(4.1.6) ^D( fe (-).fe(0) = | fo(|x|- 1 *p)c/x+M^)-^D(p,p) 

and therefore 

(4.1.7) H w > ^ H M/£ (x,) + /< A/ (x)-^D(p,p)-C5- 1 A/ 

1</<A/ 

on (Jf 2 (R 3 ,C q )) N where l/l/ £ is the "smeared" potential: 

(4.1.8) W e {x) = V(x) -(p £ * |x| _1 * p. 

Note, that the smeared potential doesn't depend on x and is defined via p 
rather than p^. 

22 ) That uniformly charged sphere S(0, r) creates potential v(x) — — q min(|x| _1 , r _1 ) 
where q is the total charge of the sphere. 
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Also let us define 

(4.1.9) A/ X (x 2 x N ) := (Xx*<Ps){xj) 

2<j<N 

and 

(4.1.10) N x := J p(y) X x(y)dy 

with xx(y) := x(x, y), x e ^°°(m 5 ). 

Furthermore, let us consider function 9 e ^^(IR 3 ) such that 

(4.1.11) O<0<1, |V Q ^| < c Q b~ H Va. 
Consider 

(4.1.12) ^ := '/ (^( x ^) -p(yWM)0(x)x(x,y)c/x</y|. 
Obviously 

J < J pi?W)|Xx * 0e " Xx(y)|0(x) c/xc/y 

+ N J |V|/(x, x 2 x/v)| 2 | A/ X (x 2 x w ) - N x \6{x) dxdx 2 ■ ■ ■ dx N 

< CNe s \\V s y +1 X \\^Q 

+ {N J |v|/(x,x 2 x N )\ 2 \N x -N x \ 2 6{x)dxdx 2 ---dx N } 

where Qy\-, •) is the quantum correlation function, 

(4.1.13) £?JW) : = A/(A/-1) y |V(x,y,x 3 x w )| 2 c/x 3 • • • dx N , 

(4.1.14) y ^\x,y)dy = (A/ - l)p„(x), 

(4.1.15) = 0^ : = J e{x) P M,(x)dx 
and we used Cauchy-Schwartz inequality. 
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Since 

(4.1.16) A/ X (x 2 x N )-N x = I ( <t>s{y-Xj)-p(y))x(x,y)dxdy, 

2<j<N 

we again from Cauchy-Schwartz inequality conclude that 

(4.1.17) \N X - N x \ 2 < C\\V yX \\%2 ■ K/v-i(x 2 x N ). 

Note that 

(4.1.18) WHxj^^^xj)^) = 

J J 9{x)p^{x)dx + J \V9 1 2\ 2 {x)p M ,{x)dx 



l<j<N 



and then (4.1.7) yields that 



(4.1.19) E/v J 9(x)pv(x) dx > - J \V9 1 2\ 2 {x)p^{x)dx 

+ £ (H w .{x k )oHxj)v.e{xi)*) 



l<j,k<N 



+ J K N {x)6{x 1 )\^{x 1 x N )\ 2 dx 1 ...dx N -^D(p,p)Q-Ce~ 1 NQ. 



Also note that 



(4.1.20) The sum of (N — 1) lowest eigenvalues of H We on H is greater than 
(WV + N^AV -*/)). 

Then the second term in the right-hand expression of (4.1.19) is bounded 
from below by (yN + Ni(/-/ M / e — u))Q, while the left-hand expression is E N Q. 
Therefore assembling terms proportional to we conclude that 

(4.1.21) se+ / \ve^\ 2 p^dx>Y,{H w ^{xj)^,e^{xj)^) 

J j 

+ N J K N { Xl x N )9{x 1 )\^{x 1 x N )\ 2 d Xl ■ ■ ■ dx N 
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with 

(4.1.22) S:=E N -uN- N l (H We - v) + ^D{p,p). 

Due to non-negativity of operator D 2 , the last term in (4.1.21) is greater 
than —CTQ with 

(4.1.23) T= sup l/l/; 

supp 6 

so we arrive to 

(4.1.24) N J M*i x w )^(x 1 )|v|/(x 1 x N )\ 2 d Xl ---dx N < 

C(S + 7 + e- 1 A/)0 + P 

with 

(4.1.25) P = / l^^l 2 ^^- 
Combining this with (4.1.13), (4.1.17) we conclude that 

(4.1.26) J< Csup||V y Xx||^( R 3 ) ((S + £- 1 A/+ 7)0 + p) 2 0^ 

+ C£/V||V y x|Uoo0 

due to obvious estimate 

(4.1.27) Ka/-i(*2 x N ) < 2K N ( Xl x^ + e^N. 

Now we want to estimate S from above and for this we need an upper 
estimate for E w . Recall that due to arguments of subsection 1.2 S < CQ 
provided we manage to prove that expressions (1.2.7)-(1.2.10) satisfy the 
same estimates as before if we plug W £ instead of 1/1/. 

So we need to calculate both semiclassical errors (which are calculated 
exactly as for 1/1/ = W TF ) and the principal parts, and in calculations of 
D(P'(W £ + p TF , P'(W £ + p TF ) and D(p £ - p TF , p £ - p TF ) an error 
is 0(e 2 Z 3 ) due to estimates 

(4.1.28) \D a (W-\N £ )\< 



€ 2 r 6 ~ H Va as £>eZ—3j^r, 

e 2 r 6 ~ H Va:|a|<| as £ x r. 
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Then one can prove easily that the sum of these two expressions does not 
exceed CQ + CZ 3 e 2 + CZ 2 e, and this estimate cannot be improved. Choosing 

2 5 

e < Z 3 we estimate these two terms by CZi. 

Under this restriction an error in the principal part of asymptotics of 
J e(x,x,\)dx does not exceed CZe 2 which is less than the semiclassical 
error. Then 5 < CQ. 

So, the following proposition is proven: 

Proposition 4.1.1. If9,x are as above then 

(4.1.29) J=\ J(/i\x,y)-p{y)p v {xfje{x)x{x,y)dxdy\< 

csu P ||v y x x ||^ 2(R 3 ) (((? + £- 1 A/+ r)ie + piei) + CeN\\v yX y^e 

with = 0x|/ defined by (4.1.15) and T, P defined by (4.1.23), (4.1.25) and 

2 

arbitrary e < Z ? . 

4.2 Excessive negative charge 

Let us select 9 = 9 b 

(4.2.1) supptf C {x : £(x) > b}. 
Note that H N W = E N W yields 

(4.2.2) E w J pv{x)e{x)8{x) dx = J2 0", e(xj)6(xj)H N V) = 

j 

(e(xj)hHxj)v,e(xj)hkxj)H N v) - \\v(oH*j)t(xi)*)n 2 

j j 
and isolating the contribution of j-th electron in j-th term we get 

(4.2.3) E N J pv{x)i{x)e dx > E N ^ J PxV (x)e(x)9(x) dx+ 

£oM*y)*foo(-v(*;)+^ 
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due to non- negativity of operator D^. 

Now let us select b to be able to calculate the magnitude of 0. Note 
that 

(4.2.4) | J e(x)(pv(x) - p(x)) dx\ < CD(p„-p,pv-p)5||V^|| ~ 

CD(pm — p, — p)^b^ < CQ^bi 



dx 



and 

(4.2.5) J 6{x)p(x) 
as long as 

(4.2.6) Z-i < b < C(Z- N)l* 

x _i 
because p x |x|~ 6 as Z '~3 < |x| < c (Z — A/) + 3 . Note that the right-hand 

expression of (4.2.5) is larger than the right-hand expression of (4.2.5) as 

b > C (? _ 7. Therefore let us pick up 

(4.2.7) b:=e Q-^; 
it does not conflict with (4.2.6) provided 

(4.2.8) N>Z-C Q$ 
and then 

(4.2.9) 9x<?'. 

Then the same arguments imply that the last term in the right-hand expres- 
sion of (4.2.3) does not exceed Gf> _1 0; using inequality 

(4.2.10) J p^,{x)i{x)d{x)dx > bQ 
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we conclude that 

(4.2.11) b\ N Q b < J 9(x)V(x)e{x) P M,(x)dx 

- J ^ 2) (x,y)£(x)|x - y|-^(x) dxdy + Cb^Q = 
= J 0{x)V{x)e{x)pv{x)dx 

- J £\x,y)£(x)\x-y\- 1 (l-e(y))e(x)dxdy 

- f g^(x,y)£(x)\x-y\- 1 9(y)9(x)dxdy+Cb- 1 Q. 

Denote by Xi, X 2 , and X3 the first, second and third terms in the right- 
hand expression of (4.2.11) respectively. Symmetrizing X 3 in the right-hand 
expression of (4.2.11) with respect to x and y 

X * = -\J Q ( i\x,y){((x) + i(y))\x-y\- 1 6(y)6(x)dxdy 

and using inequality £(x) + £(y) > min 7 (|x — xj\ + \y — Xj|) > |x — y| we 
conclude that this term is less than 

(4.2.12) - l -f £\x,y)6(y)e(x) dxdy = 

-l(N-l) J pxv (x)9(x)dx+ 1 - J Q ( £\x,y)(l-6(y))6(x)dxdy. 

Here the first term is exactly — |(/V — 1)0; replacing Q^\x,y) by p(y)p\v(x) 

we get 

(4.2.13) \ J (l-6(y))p(y)dyxQ 
with an error 

(4.2.14) \ J (^ 2) (x,y) - p(y)p^x)) (1 - 0{y))0(x) dxdy 

which we estimate using proposition 4.1.1 with x( x >y) = 1 — @(y)- Then 
||V y x x ||jf2 x />§, ||V y x||^foo x b-\ T x b~\ P x b^Q and picking up 
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e = Z 3 we conclude that an error (4.2.14) is less than Cb^Q^Q x CQ^Q 
and we conclude that 

(4.2.15) X 3 < ~{N-Z)+e+ CQ 3 iQ 

because J p9 dy x (?t and J p(y) dy = min(Z, A/). 
On the other hand, 

(4.2.16) X 2 <- J Q <t\x,y)e(x)\x - yl-^l - 9(y))9(x)dxdy 
with 9 = 9t,(i- e ) and replacing g\, (x,y) by p(y)pvj/(x) we get 

(4.2.17) -y"^)A<y^M|x-y|- 1 (l-fl(y))fl(x)dwfy 



, 3 



and we estimate an error in the same way by CQiQ. 
Therefore 

Xi +X 2 < y 0(x)\/(x)£(x)pv(x) cbc- 

y M*)t{x)p(y)\x - y\'\ l - o(y))o{*) dxd v + cq 3 q = 

J 9{x)W{x)e{x) Pxv {x)dx 

+ J p{y)pM,{x)£{x)\x - y\- 1 9{y)9{x) dxdy + CQ^Q 

due to V — W = |x| _1 * p; since W£ < Cb~ 3 we can skip the first term in 
the right-hand expression. Furthermore, as 



/ 



P(y)k-yr%)c/yx0x^ 



we can skip the second term as well. 

Adding X 3 and multiplying by _1 we arrive to 

(4.2.18) b\ N <~(N-Z)+ + CQ 3 7 

which implies immediately 
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Theorem 4.2.1. Let condition (2.2.15) be fulfilled, 
(i) In the fixed nuclei model let l/v > 0. Then 

3 5 f 1 as a < Z^, 

(4.2.19) (N-Z) + <CQy = CZU , ~ , 

(wj In particular for a single atom and for molecule with a > Z~i +Sl 

(4.2.20) (A/-Z)+<zH ; 

(Hi) In the free nuclei model let l N > 0. Then estimate (4.2.20) holds. 
4.2.1 Estimate for ionization energy 

3 

Recall that as N < Z we assumed that N > Z — CQ? (see (4.2.8)) and 
b = Q~i . Then (4.2.18) also implies In < CQ? and we arrive to 

Theorem 4.2.2. Let condition (2.2.15) be fulfilled and let N > Z — C Zh 

Then 

(i) In the framework of fixed nuclei model 

20 

(4.2.21) l w <CZs. 

(ii) In the framework of free nuclei model 

(4.2.22) Xv < Z 2 i- 5 '. 

Remark 4-2.3. (i) Classical theorem of G. Zhislin [Zh] implies that the 
system can bind at least Z electrons; the proof is based on the demonstration 
that the energy of the system with N < Z electrons plus one electron on the 
distance r is increasing as r — > +oo because potential created by the system 
with N < Z electrons behaves as (Z — N^x]' 1 as |x| — > oo; 

(ii) In the proof of theorem 4.2.1 (iii) and 4.2.2(h) we note that tearing off 
one electron in free nuclei model is easier than in the fixed nuclei model. 

The following problem looks relatively easy: 
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Problem 4.2.4. In statements (i), (ii) of theorems 4-2.1 and 4-2.2 get rid 
off condition (2.2.15). 



5 Positively charged systems 

In this section we consider the case of positively charged system with 

3 

Z — N > C (?7 with sufficiently large C . 

First let us find asymptotics of the ionization energy; the principal term 
will be — v but we need to estimate a remainder. 



5.1 Estimate from above for ionization 
energy 

As M = 1 construction is well-known: let 
us pick up function 9 such that 9 = 1 as 
\ x ~ y m | > ~f — b and # = as |x — y m | <? — 2b 
where r is an exact radius of support p TF . 
Here fccf. As M > 1 let us pick instead 

(5.1.1) 9(x) = f 2 (v- 1 [W(x) 



']) 




W(x) 



1-0 = 1- 

supp 9 



where f G ^°°(R), supported in (— oo, 2) and 
equal 1 in (— oo, 1), v < v. 
We will assume that 

(5.1.2) a> rx G(Z - 

with sufficiently large C\\ we will discuss dropping this assumption later. 
Then as we know that p TF is supported in cr-vicinity of nuclei, we conclude 
that "atoms" are rather disjoint. 

One can see easily that then as W + v w 0, 

s-5. 



(5.1.3) 



IVH/I 



then the width of the zone {x : 0< W(x)+v < 2v} is x t;|VH/| _1 
and 



b = vr 5 



(5.1.4) 







TF 



:= J 0(x)p TF 



dx x v 2 x br = v 2 r 
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while ||V0|| x b = v *r 2 and therefore to ensure that has the same 
magnitude (5.1.4) we pick up the smallest v such that v^T 1 > Cv~^?~^Q^ 
i.e. 



(5.1.5) v:=C 2 T-^Q^ ( b= QffQs)- 
then 

(5.1.6) 9xr^(Js 



Then (4.1.15) is fulfilled. Note that v < u = T~ 4 iff (Z - N) > Qi exactly 
as we assumed. 

Then (4.2.11) is replaced by 

(5.1.7) \ N J e(x)pv(x)8(x)dx< J 9{x)V(x)e(x) P M,{x)dx 

- J Q^\x,y)£(x)\x - y|-^(x) dxdy + Cb~ 2 rQ 

where Cb~ 2 ?Q estimates the last term in the right-hand expression of (4.2.3). 
Then 

(5.1.8) \ N J e{x)pv(x)6{x) dx < J 9{x)V(x)e{x)pv(x)dx 

- J (ti\x,y)-p*(x)p(y)y(x)\ X -y\- 1 6(x)dxdy 

- J P*(x)p(yy(x)\x - y\~ X 0{x) dxdy + Cb- 2 ?Q. 
Let us estimate from above 

(5.1.9) - J(Q {2 \x,y)-pv(x)p(y)y(x)\x-y\- 1 9(x)dxdy< 

- J r) - M*Mr)) (1 - w(x, y))t{x)\x - y|-^(x) dxdy 

+ J py(x)p(y)uj(x,y)e(x)\x - y\' 1 9(x)dxdy 

with u = w 7 : w = as |x — y| > 27 and u> = 1 as |x — y| < 7, 7 > b. 

To estimate the first term in the right-hand expression one can apply 
proposition 4.1.1. In this case ||V y x||j&? 2 x C?^~2, || V y x||j^°° x ?7~ 2 and 
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plugging P = b 2 and T = v, e = Z 3 we conclude that this term does 
not exceed 

(5.1.10) Cr( 7 -^ + CZi 7 " 2 )0, Q = max((J, Z§). 
Consider the second term in the right-hand expression of (5.1.9). Note that 

J p(y)u(x,y)\x - y| _1 dy < C{?'^^ + v^ 2 ) 

since p(y) x (H/(y) + and |Vl/l/| x r~ 5 ; then this term does not exceed 
C(7~t 7 5 +-u§7 2 )r0. 

Adding to (5.1.10) we get 

/ _1 -il 2 _ 7 15 7 3 <j\__ 

CI7 2Q2 + CZ37 + r 2 72 -|- W 27 z 1/-0. 

— 1 3 3 — 2 

Optimizing with respect to 7 = Q$v~$ we get OuioQsrO and (5.1.8) 
becomes 

(5.1.11) (l/v + v) J e{x)pv(x)9{x) dx < 

J 6(x)(W(x) + v)i{x)pv{x) dx + CvTo Ql?Q < 

C(v + v^o Q^)?Q. 

where we took into account that V — \x\~ l * p = W, and that W + vv on 
supp6>, 

Since a factor at (l/v+f) hi the left-hand expression of (5.1.11) is obviously 

3—2 

x 70 we arrive to (l/v + z/) < Cv + fie Q5. 

Recalling definition (5.1.4) of v we arrive to an upper estimate in Theo- 
rem 5.2.1 below: 

(5.1.12) \ N + v < Cv = C(?5(Z- N) n. 

Really, this is true if (? > Z? (we are interested in the general approach) 
and in this case v > fioQs. On the other hand, if (? = Z^aZa) -5 , 3Z3 > 1, 

_ 1 12 

then we get an extra term C(?5o(Z — N)^Z^ but we can skip it decreasing 
unspecified exponent 5 > in the definition of Q. 
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Remark 5.1.1. (i) We will prove the same estimate from below in the next 
subsection. 

(ii) Note that the relative error in the estimates is v/v = (Z — N)~iQ*. 

(iii) In the proof we used not assumption (5.1.2) itself but its corollary 
(5.1.3). If we do not have such assumption instead of equality (5.1.4) we 
have inequality TF >V2? 7 (but probably equality still holds). However we 
have a problem to estimate |0 — TF |: namely we need estimate the first 
factor in the product || W 5 ||D(p — p TF , p — p TF ) 5 . 

Let us select f in (5.1.1) such that \ f'\ < cf 1 ^ 5 ^ 2 with arbitrarily small 5 > 0. 
Then 

b 2 ||V^|| 2 <C I 6 1 - & dx<C{j 9dxy- s C{[ ldx) 5 < 

*/ Z o Z iJ Z 

C(0 TF ) 1 -^-i( 1 - 5 )r 35 

with Z = supp W6 and b = ?^Q^> and an error is less than e0 TF provided 

v = C 2 ?-TQl x {? 7 Q)- 51 
which leads to a marginally larger error. 

Estimate P could be done in the same manner but here slight increase of it 
does not matter. 

(iv) The same arguments of (iii) could be applied to the proof of the lower 
estimate in the next subsection despite rather different definition of by 
(5.2.6); 

(v) In view of Theorem 5.3.1 stable molecules do not exist in the free nuclei 
model as Z — N > CQi and in atomic case I a/ = l/v- 

Problem 5.1.2. Consider f such that \ f'\ < fg(f) where g~ 1 (t) G Jz? 1 and 
further improve remainder estimate without assumption (5.1.3). 
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5.2 Estimate from below for ionization 
energy 

Now let us prove estimate I a/ + v from below. Let W — W/v(xi, ... , x/v) be 
the ground state for N electrons, ||W|| = 1; consider an antisymmetric test 
function 

(5.2.1) vj> = vj>(xx x/y+i) = v|/( Xl x N )u{x n+1 )- 

^ ^(Xi Xy_i, X/y+i , X; + i X N )u(Xj) 

l<j<N 

Then 

IM| 2 < (Ha/+i^>) = N(H N+1 Vu,V) = 
N(H N Vu,V) + N(H v , XN+1 Vu,V) + A/( ^ |x / -x /v+ i|- 1 M/£v,i) = 

1</<W 

(E w - ^llVf + N(H w+u , iXN+1 Vu, vj>) 

+ A/<( ^ |x / -x A/+1 |- 1 -(\/- H/)(x w+1 ))vl/ty>) 

1<;'</V 

and therefore 

(5.2.2) A/ _1 (l w +i + i/OII^II 2 > (H w+ ^, X/v+1 V|/u, V) 

1</<A/ 

with v' > v to be chosen later. One can see easily that 

(5.2.3) AT 1 !!^!! 2 = ||V|| 2 • ||u|| 2 - 

A/ J V|/(xi x A/ _i,x)v|/ t (xi x /v _i,x)iy(y)iy t (y) cbci • • • c/xoy 

where t means a complex or Hermitian conjugation. 

Note that every term in the right-hand expression in (5.2.2) is the sum 
of two terms: one with W replaced by ^(xi, ... , xn)u(xn+i) and another with 
M> replaced by — A/v|/(xi, ... , x/v_i, x N+1 )u(x N ). We call these terms direct and 
indirect respectively. 

Obviously, in the direct and indirect terms u appears as |tv(x)| 2 dx and 
as u{x)u\y) dxdy respectively multiplied by some kernels. 
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Recall that u is an arbitrary function. 
Let us take u(x) = 82(x)<f)j(x) where <f)j 
are orthonormal eigenfunctions of H w+U and 
6>(x) is b-admissible function supported in 
{x : — v > W(x) + v > |z/} and equal 1 in 
{x : | — 2v > W(x) + v > |z/}, satisfying 
(4.1.11), and b = vT 5 . 

Let us substitute it into (5.2.2), multiply 
by ^(XjL^ 1 ) and take sum with respect to j 
we get the same expressions with |u(x)| 2 dx 
and u(x)u\y) dxdy replaced by F(x,x)dx 
and F(x, y) dxdy respectively with 



W(x) 



v 



V 



1-0 = 1- 

supp0 



(5.2.4) 



F(x,y) = J V 9(A/.- 1 )c/ A e(x,y,A). 



Here (p(r) is a fixed < ^' 00 non-negative function equal to 1 as r < \ and equal 
to as r > 1 and L = v' — v = 6v. 

Under described construction and procedures the direct term generated 
by A/- 1 !!^!! 2 is 



(5.2.5) 



J 9{x)cp{\L~ 1 )d x e{x,x,\)dx 



and applying semiclassical approximation we get 

(5.2.6) 0v := J ifiiXr 1 ) d x P' B {W + is - X) dx 

and under assumptions (2.2.15) and (5.1.2) 23 ) the remainder estimate is 
Ch~ l ? 2 b~ 2 = Cv2? 2 b~ 1 = Cv~2?~ 3 ; one can prove it easily by partition 
of unity on supp$ and applying semiclassical asymptotics with effective 
semiclassical parameter h = l/(v2b) = v~2?~ b . 

On the other hand, indirect term generated by A/ -1 !!^ 2 is 

(5.2.7) -N J ^(x)^( y )v|/( Xl x A/ _ 1 ,x)M/t(x 1 x /v _ 1 ,y)x 

F(x, y) dxdydxi ■ ■ ■ dx^j-i 



23 ) Or rather its corollary (5.1.3). 
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and since the operator norm of F(., ., .) is 1 the absolute value of this term 
doesn't exceed 

(5.2.8) N J 0(x)|V(xi x N _ 1 ,x)\ 2 dx = J 0(x)p„(x) dx < 

J 9(x) P (x)dx+ cQ^wve'm 

where p TF = on supp# and under assumption (5.1.2) 23 ) ||V0^|| x b~^T x 

>— i — z - 
C 2 r 2 - 

Recall that P'(W TF + v) = p TF . We will take v' = v + L large enough to 
keep 0vi/ larger than all the remainders including those due to replacement 
W by W Tf and p by p TF in the expression above. One can note easily that 

(5.2.9) Gxu x h~ 3 x b~ 2 r 2 x t^br 2 x vh 7 . 
Therefore 

17 1 5 

(5.2.10) As v = Qr^Qe and Z — N > CqZ? the total expression generated 
by A/ -1 !!^!! 2 is greater than e0 with = v*? 1 . 

Now let us consider direct terms in the right-hand expression of (5.2.2). 
The first of them is 

(5.2.11) -J 9 1 i(xM\L- 1 )d x (H w+u ,J 1 i(x)e(x,y,\)) y=x dx = 
- J 0(x)^(A/.- 1 )of A (H M/+ ^, x e(x,y,A)) y=x ofx 

- l -J ^\L- l )[[H w ,d^},d^]d x e{x,x,X)> 

J 6{x)(v' -v -\) V {\L- 1 )d x e{x,x,\)dx- C J | | 2 e(x, x, v')dx. 

Note that the absolute value of last term in the right-hand expression of 

(5.2.11) does not exceed Or 1 ? 2 /. 2 x Cv2?~ 3 <c CvQ. 
The second direct term in the right-hand expression is 

(5.2.12) -J 8{x)(pv* Ix]- 1 - {V- W){xfjF{x,x)dx = 

-D(p w -p,0{x)F{x,x)) > 
- CD(p M/ -p,p^-p)5 ■ D(^F(x,x),^F(x,x))) 2 > -CQh^vi 
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provided V — W = |x| _1 * p with D(p — p TF , p — p TF ) < CQ; the absolute 
value of this term is <C vQ. 

Further, the first indirect term in the right-hand expression of (5.2.2) is 

(5.2.13) -N J 02(y)V(xi x w _i, x)M/ t (x 1 x /v _ 1 ,y)x 

V?(A/. _1 ) c/ A (H M/+iy /, x 6'5(x)e(x,y, A)) dxdydx x ■ ■ ■ c6c w _i = 

-A/y 05(y)05(x)V(x! x A/ _ 1 ,x)M/ t (x 1 x w _ 1 ,y)x 

<p(A/. _1 )(z/' — z/ — A) d\e(x, y, A) dxdydx! ■ ■ ■ dx^-i 

-N J 05(y)V(xx x A/ _ 1 ,x)M/t( Xl x /v _ 1 ,y)x 

(/?(A/. _1 )[H^ X , 0^(x)] c/ A e(x,y, A) dxdydx 1 ■ ■ ■ dx N _ x . 

Note that one can rewrite the sum of the first terms in the right-hand 
expressions in (5.2.11) and (5.2.13) as <p(Ay/. -1 )(i/ — z/ — Aj)||Vy|| 2 with 

V/(xi, ... ,x w _i) := J ^(xx, ... ,x A/ _i,x)6»5(x)0 i (x) c/x 

and therefore this sum is non-negative. 

One can prove easily that the absolute value of the second term in (5.2.13) 
is less than 

Cvh' 1 J pv{y)6*(y) dy < Cv-^?- 5 Q < t>0. 

Therefore 

(5.2.14) The sum of the first direct and indirect terms in the right-hand 
expression of (5.2.2) is greater than — CvQ. 

Finally, we need to consider the second indirect term generated by the 
right-hand expression of (5.2.2) 

(5.2.15) - f(J2 |y-*l -1 -(v-M0(y))x 

J l<i<N 

v|/(xi, ... ,x A /)v|/ t (xi, ...,x N - lt y)e*(x N )6*(y)F(x N ,y) dx x ■ ■ ■ dx N dy = 
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0Hx N )9 1 2(y)F(x N , y) dx x ■ ■ ■ dx N dy 

^2 \y- x i\~ 1 -\y\~ 1 *Px(y) s )v(xi, ...,x N )¥(x 1 , ... ,x w _i,y)x 



Ki<N 



9~ 2 {x N )92(y)F(x N ,y) dx l --- dx N dy; 

recall that is a smeared density, x = (x 1: ... ,x N ). 

Since {y]' 1 * Px{y) — (V — W)(y) = | | 1 * (p K — p), the first term in the 
right-hand expression is equal to 

(5.2.16) J ^(x A/ )v|/(x 1 x w )x 

D y(px(y) - p(y), F(x N ,y, A)05(y)V(xi, ... .xw-i.y)) d*i • • • c6c w ; 
and its absolute value does not exceed 
(5.2.17) 

(/V J D(p,(-)-p(-),Px(-)-p(-))l v l ; (^i x^l^x^t/xi-.-t/xA/yx 

AH |^D y (f(x n , y, A)0*(y)V(xi x N . lt y), 

F(x N ,y,X)9 1 2(y)^f(x 1 x/v_i,y)) dx 1 ---dx N ^j . 

Due to estimate (4.1.24) and definition (4.1.4) as the first factor in 

i 

(5.2.17) does not exceed (((? + T + e~ 1 N)Q + P) 2 where we assume that 
e < Z~i and x b~iQi? is now an upper estimate for J 0{y)p^{y) dy-like 
expressions; due to our choice of v it coincides with = v^l' ' . 

Then according to (4.1.25) P x Cb~ 2 Q < QQ and according to (4.1.23) 
T <^ Q and therefore in all such inequalities we may skip P and T terms; 
so we get C(Q + e 1 A/)l©l . 

Meanwhile the second factor in (5.2.17) (without square root) is equal to 

AT 1 / /--V(A/.- 1 y(A / /-- 1 )|y-z|- 1 e(x A/ ,y,A) ^(y)v|/(x 1 x A/ _ 1 ,y)x 

e(x N , z, A') e 1 2(z)^(x 1 x w _i, z) dydz dx 1 --- </x w _i c/x w c/Ac/A'; 
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after integration to x/v we get instead of marked terms e(y, z, A) (recall that 
e(., ., .) is the Schwartz kernel of projector and we keep A < A') and then 
integrating with respect to A' we arrive to 



N 



'J \y-z\~ 1 F{y,z)9 1 2{y)^{x 1 x /v _ 1 ,y)x 

6'5(z)M/ t (xi x/v_i, z) dydz dx x ■ ■ ■ dx N _i\ 



where now F is defined by (5.2.4) albeit with (p 2 instead of ip. This latter 
expression does not exceed 

(5.2.18) N 1 JJ |y-z|- 1 |F(y,z)|^(y)|M/(x 1 x N _ lt y)\ 2 x 

dydz dxi ■ ■ ■ c/x/v-i. 

Then due to proposition 6.C.1 f \y — z\~ 1 \F(y, z)| dz does not exceed 
Cb^ 1 ^ 1 x vz, and thus expression (5.2.18) does not exceed CZ~ 2 vzO 
and therefore the second factor in (5.2.17) does not exceed CN~ l v^Q^ and 
the whole expression (5.2.17) does not exceed 

C{Q + e- 1 N) 1 2Q 1 2 x ArVe? = CN~\Q + e~ x N)W*Q 

and then 

(5.2.19) As e > Z~ l v~^ the first term in the right-hand expression of (5.2.15) 
does not exceed CvQ. 

Further, we need to estimate the second term in the right-hand expression 
of (5.2.15). It can be rewritten in the form 

(5.2.20) I U(xi.yW><i xn)^(xi x /v _ 1 ,y)^(x A/ )^(y)x 

l<i<N J 

F(x N ,y) dx x ■ ■ ■ dx N dy 

where l/(x,-,y) is the difference between potential generated by the charge 
5(x — X,) and the same charge smeared; note that U(x,,y) is supported in 
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{(x,-, y) : \x; — y\ < e}. Let us estimate the ;-th term in this sum with /' < N 
first; multiplied by N(N — 1), it doesn't exceed 

(5.2.21) 

N^j \U{x„y)\ 2 \Ai{x l x N )\ 2 9 1 Hx N )9 1 i(y)\F(x N ,y)\dx 1 ---dx N d y yx 

N[ [ cu{x i ,y)\^{x 1 x N - 1 ,y)\ 2 6i{x N )6s{y)\F{x N ,y)\dx 1 ---dx N dy 



here u> is e-admissible and supported in {(x,-,y) : |x,- — y| < 2e} function. 
Due to proposition 6.C.1 in the second factor J 62 (x N )\F(x N , y)| dx N < C 
and therefore the whole second factor doesn't exceed 

(5.2.22) C^j 62-(x)oj(x,y)Q { £\x,y)dxd y y 

where we replaced x,- by x. According to Proposition 4.1.1 in the selected 
expression one can replace g\, (x,y) by pvj/(x)p(y) with an error which does 
not exceed 

c(sup||V y Xx||^ ( R3 ) ((? + £- 1 A/)^ + CeN\\V yX y°°)e 

which as we plug sup x || V y Xx||_sf 2 (R 3 ) x W^yX\\jf°° x e~ l becomes CNQ. 
Meanwhile, consider 

(5.2.23) J \U(x h y)\ 2 e-2(y)\F(x N ,y)\dy. 
Again due to proposition 6.C.1 it does not exceed 

Cv\ J \U(x i ,y)\ 2 9 1 2(y)(\x N -y\v 1 2+iy s dy 

and this integral should be taken over 6(x;,e), with \ U(xj,y)\ < |x,- — y| , 
so (5.2.23) does not exceed Cev2u'(xj, x^) with u/(x,y) = (l + V2~\x — y|) s 

(provided e < which will be the case). Therefore the first factor in 
(5.2.21) doesn't exceed 

(5.2.24) Ce^v^J 9*{x)u' (x, y)^\x, y) dxdy) 2 . 
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Therefore in selected expression one can replace Q^\x,y) by Pw(x)p(y) with 
an error which does not exceed what we got before but with e replaced by 
v~2 i.e. also CNQ. 

However in both selected expressions replacing (x, y) by p\j/(x)p(y) we 

1 3 

get just 0. Therefore expression (5.2.21) does not exceed Ce^v^ZQ which 
does not exceed CvQ provided e < Cv^Z^ 2 . 

2 

So, we have two restriction to e from above: the last one and e < Z~i 
and one can see easily that both of them are compatible with with restriction 
to e in (5.2.19). 

Finally, consider term in (5.2.20) with /' = N (multiplied by A/); 

(5.2.25) nJ U{x N ,y)M Xl x N )\ 2 e^{x N )e^{y)F{x N ,y)dx 1 - ■ ■ dx N dy 

due to Cauchy inequality it does not exceed 

(5.2.26) N^J |x A/ -y|- 2 |M/(x 1 x A/ )| 2 ^^x A/ )^5(y) d Xl ■ ■ ■ dx N dy^j 2 x 

N(J \F(x N ,y)\ 2 \V( Xl x N )\ 2 e*{x N )6*(y) d Xl ---dx N d y y 

where both integrals are taken over {|x/y — y| < e} and integrating with 
respect to y there we get that it 

Ce^Qi x v^e^Q^ = Cv*e 2 Q vQ. 

Therefore the right-hand expression in (5.2.2) is > — CvQ and recalling that 
v 1 — v = 0(v) we recover an lower estimate in Theorem 5.2.1 below: 



(5.2.27) \ N + is > -Cv = C(?s(Z- N) 



Combining with a lower estimate (5.1.12) and recalling estimate (3.1.5) 
for Q we arrive to 

Theorem 5.2.1. Let conditions (2.2.15) and (5.1.2) be fulfilled and let 
N <Z - C Q 3 i with Q < CiZl . Then 

In the framework of fixed nuclei model under assumption (5.1.2) 



(5.2.28) + H < C(Z- N) is Z 



17 

' 18 



1 as a < Z 3, 

AT 5 + (aZ5)~ 5 as a > Z '"i . 
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5.3 Estimate for excessive positive charge 

To estimate excessive positive charge when molecules can still exist in free nu- 
clei model we apply arguments of section 5 of B. Ruskai and J. P. Solovej [RS]. 
In view of Theorem 3.4.2 it is sufficient to consider the case 

(5.3.1) a = min \xj — x k \ > CqT. 

j<k 

Therefore in Thomas- Fermi theory p TF is supported in separate "atoms" . 

Consider a-admissible functions m (x), supported in B{y m , ^a) as k = 
1, ... , M and in {|x — y m | > |a Vm = 1, ... , M} as k = 0, such that 

(5.3.2) 6l + ... + 9 2 M = 1. 
Then for the ground state U' 

(5.3.3) E N = (HV|/, (0 a H0 a V, v|/) - ^ II V^^|| 2 

a a,k 

with the sum over of (M + l)-cluster decompositions a = (cyq, ... , a/w) of 
{1, ... , A/} and 6 a {x) = Uo<k<M Hiea k °k{xi)- Then for any given a 

(5.3.4) H= H - k + J a 

0<k<M 

with the cluster Hamiltonians H ak involving only potential of /c-th nucleus 
(no nucleus potential as k = 0) and only electrons belonging to a k and such 
that 

(5.3.5) H ak > E at (N k (a),Z k ), H ao > 
and with the intercluster Hamiltonian 

(5.3.6) J a = ^-^k-y^r 1 

l<k<M i(£a k 

+ Y1 E \ Xi ~ x j\~ l + ^2 z k z i\yi -y/cT 1 - 

k<l iea k ,jeai k<l 

Let us note that 

(5.3.7) £^4 = E J > 



kl 

0<k<KM 
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with J k i given by (32)-(33) of Ruskai-Solovej [RS] as k, I > and k = 
respectively: 



(5.3.8) J k i = Z k Z,\y k -y/]- 1 - Z, ^/(x,) 2 ^- 



y* 



-1 



Z^^x,) 2 ^ - y/r 1 + Y,d k {x i fe l {x j f\x i -x ] \-\ 

and 

(5.3.9) j / = ^ ^( x <) 2 (- z /l x / - yl -1 + E 0'(*/) 2 l*«- - ■ 

Then we recover (35) Ruskai-Solovej [RS] 

(5.3.10) (J W M/, VI/) = Z fc Z,|y* - y^ 1 - Z, J p^{x)9 k {xf\x - y^ 1 dx- 

Z k J pv{x)6 l {x) 2 \x-y k \- 1 dx + J g i2 \x,y)9 k {x)9 l (y)dxdy. 

Applying proposition 4.1.1 and estimate (3.4.1) (replacing first 9 k with 
k = 1, ... , M by 9 k supported in B(x k , c?) and estimating an error), we 
conclude that 

(5.3.11) J pv{x)9 k {x) 2 \x - y,!" 1 dx = 

( J p TF (x)9 k (x) 2 \x - y,!- 1 dx + O(V)) \y k - y^ 1 , 



(5.3.12) J P v{x)9 k {x) 2 dx = Nj f + 0{Y), 
with 

(5.3.13) Nj F = J p TF {x)9 k {x) 2 dx, Y :=Qh^ 
(compare with (36)-(37) of Ruskai-Solovej [RS]) which yields 

(5.3.14) f Pxv (x)(l- k {x))dx<CY 



Kk<M 
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and we prove that (5.3.12) holds for 9 k as well (of Ruskai-Solovej [RS]). 

The last term in (5.3.11) is estimated by proposition 4.1.1 and estimate 
(3.4.1) and the same replacement trick: 

(5.3.15) J Q < i\x,y)e k {x) 2 e l {y) 2 \x-y\- 1 dxdy> 

J e$ ) {x,y)O k {x) 2 O l (y) 2 \x-y\- 1 dxdy > 
J p TF (x)p^y)9 k (x) 2 9 l (y) 2 \x - y\~ l dxdy- 

c(q* J pv(x)6 l (x) 2 dx+ Va-^ly.-y/r 1 
and repeating the same trick we get that it is larger than 

(5.3.16) J p TF {x)p TF {y)\x - y^dxdy - C{Z - N)Ya 1 - CY 2 



' 2 a~\ 



Then we conclude that 

(5.3.17) 
with 



(JkiV.V) > J T kl F - CNTa- 1 



(5.3.18) J J J = f p TF (x)p TF (y)9 k (x)9 l (y)\x-y\- 1 dxdy 

-Z k J p TF (x)9 l (y)\x-y k \' 1 dx-Z l J p TF (x)9 k (y)\x - y^ 1 dx 

+ Z k Zi\x k -x/| _1 

and 

(5.3.19) |(Jo/^,^)| < C(Z - A/)Va _1 + CY 2 a~ 1 
(compare with (39)-(40) of Ruskai-Solovej [RS]) provided 



(5.3.20) | J p TF (y)\x - y^O^x) dy - N T k F \x - y^ 1 



< 



C{Z- N^x-x.l' 1 
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as |x — y^| > Cr. 

Let us note that the absolute value of the last term in the right-hand 
expression of (5.3.3) doesn't exceed Ca~ 2 Y due to (5.3.12). Now stability 
condition yields 

(5.3.21) J = ^lki< CYa 2 + C(Z - A/)^ 1 + CY 2 a~ 1 . 

k<l 

This inequality, (5.3.1) and proposition 5.3.3 below yield that Z — N < 
CY = C?i<?i. Since 7 x (Z - A/)" I we arrive to (Z - N) < CQi: 

Theorem 5.3.1. Let condition (2.2.15) be fulfilled. Then in the framework 
of free nuclei model with M > 2 the stable molecule does not exist unless 

(5.3.22) Z-N<Z^- 8 . 

Remark 5.3.2. Unfortunately, we don't prove that molecules exist. We are 
not aware of any rigorous result of this type in frames of our models. 

Proposition 5.3.3. Let (5.3.1) be fulfilled. Then inequality (5.3.20) holds 
and 

(5.3.23) J >e{Z - N) 2 3 -\ 
Proof. Proof Note first that 

(5.3.24) £ TF < S(p) = J>(pJ F ) + J TF (p) < 

j 

^(pD + CiZ-Nfa-i 
j 

with p = J2j pJ F while 

(5.3.25) E^7 F )<E^(V TF ); 

j j 

then 

(5.3.26) J TF {P TF ) < C{Z - Nfa- 1 
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and using (5.3.22) we conclude that 

(5.3.27) D(p TF - p, p TF -p)< C{Z - N) 2 a-\ 
Further, 

(5.3.28) A :=J2 D(p T ^i - pj" , p TF 0j ~ pJ F ) < 

j 

D(p TF - p, p Tf -p) + CTa- 1 A 
and combining with (5.3.26) we conclude that 

(5.3.29) A < C(Z - N) 2 a~ l 

due to (5.3.1). Combining with pJ F * | ~ 1 = Nj F \x — xj| _1 for |x — xj\ > r$ 
we arrive to (5.3.20). Further, 

(5.3.30) J TF (P TF ) > J TF {p) ~ CAh^{Z - N)a- 3 2 - CA? 1 

which together with (5.3.28) and (5.3.29) yields (5.3.23). □ 

6 Appendices 

6. A Electrostatic inequalities 

We know already that there are two sources of errors in the lower estimate: 
due to electrostatic inequality (1.1.1) and semiclassical errors. For the first 
error in the case B = const E. Lieb, J. P. Solovej and J. Yngvarsson [LSY2] 
provide the (almost) perfect estimate; the reader can find the proof based 
on magnetic Lieb-Thirring inequality (and this inequality as well) in that 
paper (p. 122) which in the case of B = becomes 

Theorem 6.A.I. For the ground state W of (0.1.1) with potential (0.1.4) 
(6.A.1) J p\dx < CZ*N*(Z + N)* 

otherwise. 

In particular for c _1 A/ < Z < cN the right-hand expression doesn't exceed 
CZl 
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On the other hand, for 6 = there is a more precise inequality due to 
V. Bach [Ba] and G. Graf and J. P. Solovej [GS]: 

Theorem 6. A. 2. For the ground state of (0.1.1) with potential (0.1.4) 
(6.A.2) (HM/,M/)> 

M A ~ u ) + uN ~ \ J J |x-yr 1 |e(x,y,z/)| 2 c/xc/y- CA/f- 5 

with some exponent 5 > 0. 

We will discuss magnetic field case in more details in the Appendix to 
the next Chapter. 

6.B Hamiltonian trajectories 

We are going to prove that for 1/1/ = W TF in M = 1 case the generic 
trajectory on the energy level v is not periodic. We use some ideas from 
V. Arnold [A], pages 37-38. Recall that in this case 1/1/ = W(r) (r =\x-x 1 \) 
and angular momentum M is a motion integral. Then any trajectory lies on 
some plane and if M — \ M\ > it lies in{0<r<7} where 1/1/ is analytic 
and 1/1/(7) = -v. 

(6.B.1) Let us assume that all the trajectories on the energy level v are 
periodic. 

Then the rotation number 

T 2 Mdr 

(6.B.2) 



* = / " 



2 ^2{W(r) + v) - M 2 r- 2 



showing the increment of the polar angle over a half-trajectory should be 
7r/c _1 with k E N and should not depend on M where r x < r 2 are roots of 
(W(r) + v) — M 2 r~ 2 = v. So, should be the same for all trajectories on 
the energy level v. One can see easily that for M — > 4> tends to those of 
the Coulomb potential. So, <t> = 2n for all trajectories on the energy level v. 

Let r be a root of F(r) = 2(W(r) +u) + r 4 (W(r)) 2 = 0. One can see 
easily that F(7) > 0, F(r) — > oo as r — > and F'(r) < because 

(6.B.3) W < 0, W" + 2rW = r l (r 2 Wj) = rA W > 0. 



65 



So, the unique root exists. Then r = r$ is a circular motion with M = 
-r$W'{r ). 

Then (V. Arnold [A], problem 2 at page 37) 

<t> -> iry/W'/{3W' + rW")- 1 \ r ^ rQ = cD 

for trajectories tending to circular. However, 3W + rW" > W due to 
(6.B.3) and then <t> > n. Contradiction to assumption (6.B.1). 

6.C Some spectral function estimates 

Proposition 6.C.I. For Schrddinger operator W G and for <p G 

< ^o 5 °([— 1, 1]) the following estimate holds for any s: 

(6.C.1) \F{x,y)\ < Ch- 3 (1 + h-^x - y\y s , 

(6.C.2) F(x,y) := J 0(A) d x e{x, y, A). 

Remark 6.C.2. Surely, \AZL~ 1 is not < ^ ?0 ° even after rescaling but the proof 
works for it after rescaling. 

Proof. Let u(x,y, t) = J e~ lh lfA d\e{x, y, A) be the Schwartz's kernel of 

Fix y. Note first that J£ 2 -norm 24 ' ) of 0(/7D t )%(t)w(x)tv(x, y , r) is less than 
Ch s as x G ^^{{—t, e\) and u G supported in {|x — y| > ei} (ei = Ce) 
due to the finite speed of propagation of singularities. 

We conclude then that Jzf 2 -norm of <p(hD t )x(t)u(x)u(x, y, t) is less than 
Ch s for oo G 'To 00 supported in {|x — y| > C. 

Then J2f 2 -norm of d l t V a <f)(hD t )x(t)uj(x)u doesn't exceed Ch s . Then due 
to imbedding inequality Jzf^-norm of <p{hD t )x{t)u{x)u doesn't exceed Ch s . 
Setting t = and using this inequality and \F(x, y)\ < Ch~ 3 (due to 
Chapter 4) of [17] we get that \F(x, y)| < Ch s for |x — y| > eo; 

Now let us consider general \x— y| = r > Ch. Rescaling x—y h-> (x— y)r _1 
we need to rescale h H- /?r _1 and rescaling above inequality and keeping in 
mind that F(x,y) is a density with respect to x we get |F(x,y)| < Ch s r~ 3 ~ s 
which is equivalent to (6.C.1)-(6.C.2). □ 

24 ) With respect to x, t here and below. 
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